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PREFACE 


This unit shows you ways in which the computer can extend your 
knowledge and understanding of introductory electricity and magnetism. 
The unit presents numerical solutions for field lines and equipotentials. 
The solutions are easy to understand and yet include solutions to many 
problems which can’t be handled by analytic methods. The emphasis in 
the unit will be on ways the computer can extend the electricity and 
magnetism you understand. The unit assumes some knowledge of intro- 
ductory e and m, specifically the fields for simple charge distributions. 


Your instructor can assign chapters and exercises out of this unit in much 
the same way as out of a textbook. Lectures on the material may be 
presented, but you should be able to understand the material without 
lectures. The harder exercises can be used as the basis of projects if your 
instructor wishes. 


Fach chapter in the unit starts with a brief discussion of the physics 
discussed in the chapter and then moves on to an explanation of the 
numerical procedure used with the computer. Exercises follow with one 
exercise completely worked out in the text so that you can see what a 
sample solution looks like. This sample exercise is sometimes a problem 
for which the answer is already known and therefore provides an extra 
check on the computer method. Some exercises are marked with asterisks. 
One asterisk means the exercise is fairly demanding while two asterisks 
indicates a very challenging problem. Numerical methods will not be 
discussed much in the text; there is a section in the appendix which 
discusses the half-step integration used. 


Throughout the booklet normalized or natural units will be used. The 
units have the advantage that the numbers being calculated stay near 1; 
therefore rarely getting too small or too large for the computer. Natural 
units are often used in e and m; the units have the effect of setting 
1/(47é€9) and fto/(47) all equal to one. 


Often graphical output is useful. If you have a plotter or a CRT terminal 
available, fine. If not, terminal plotting (using the teletype to plot a graph) 
can be helpful. The appendix gives a program to convert printing programs 
{programs which type out lists of numbers) to programs which plot on the 
terminal. This terminal plotting will be used several times in the unit to 
give you several examples. 


John Merrill 
Tallahassee, Florida 


This unit was written while the author was on the faculty of the Depart- 
ment of Physics and Astronomy at Dartmouth College, Hanover, New 
Hampshire. John Merrill has used computers to support his classwork and 
research since joining the Dartmouth faculty in 1966. During the period of 
1970-72 he was part of a team of professors who developed computer uses 
for the classroom in engineering, mathematics and physics under the 
auspices of NSF funded project COEXIST. 


John has written many journal acticles on specific uses of computers in 
physics teaching, particularly for the two main physics teaching journals, 
The Physics Teacher and The American Journal of Physics. He has written 
several booklets on computer use in physics teaching including another 
unit in the Hewlett-Packard series, Quantum Mechanics. John is also the 
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author of the book, Computers in Physics, soon to be published by 
Houghton-Miffin Company. He is now the Director of the Center for 
Educational Design at Florida State University. 


Special credits go to Gregory Hughes, a recent PhD in physics from 
Dartmouth, who helped devejfop the exercises and Christine Doerr, who 
copyedited the manuscript. 
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CHAPTER ONE: ELECTROSTATICS 


INTRODUCTION 


Electrostatics introduces the basic concept of fields, a concept used widely 
in electricity and magnetism (e & m) as well as in other areas of physics. 
The purpose of this chapter is to deepen your understanding of fields by 
using the computer to map electrostatic fields in various ways. You will 
map electrostatic patterns both by using the electrostatic potential and by 
using the electrostatic field itself. 


THE ELECTROSTATIC FIELD, E(r) 


A field is a way to visualize a vector whose magnitude and direction vary 

as you move around space. Consider the Coulomb force, F,,, on a charge, 
; > 

q2, due to a charge, q,, a distance, r, away, 


> 
-> qiq r 
Fo, = ko — (=) (1) 
r r 


where r = ir| and kg defines the units. kg = 1/4m€ 9 = 9x10? in mks units; 
ky = 1 in Gaussian (cgs) units. This force means that q, reaches out over 





the distance, r, in order to create the force felt by q,. Originally, the 
electrostatic field, E was introduced into physics so that q, would be 
interacting with something (namely, the E field) right where q. was. The 
idea is that a1 produces an E field everywhere, and that q, interacts with 
the E field at qy‘s position. 


The E field turns out to be avery convenient way to handle electrostatics 
problems because you can separate the effects of q,, the so-called source 
charge, from the effects of the field upon q,. You can separate each 
problem into a part dealing with what field is produced by the source 
charges and another separate part dealing with the behavior of other 
charges placed in that field. 


The electrostatic field at some point, r= (x,y,z), is defined as the force per 
unit charge on a test charge, q., placed at r. 


Etr) = (2) 


>>. _.> 
where F(r) is the force felt by q. at the point r. 


_> 
Equations for the electrostatic field, E, for simple kinds of source charges 
are derived in regular textbooks. A short table of the fields looks like the 


following: 
> 
Charge Distribution E Field 
k, a/r? (F/r) 








Single Point Charge = the line charge density 






Single Line Charge ky Afr (r/r) 
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The total electrostatic field due to some set of source charges is just the 
(vector) sum of the fields due to each source charge individually. 


Practice Exercise: Show that the force experienced by a test charge, q2, 


placed at the point r is the Coulomb force when the source charge is a 
point charge. 


THE ELECTROSTATIC POTENTIAL, V(r) 


Another useful concept in electrostatics is the electrostatic potential, V(F). 
You can define or introduce the potential in several different ways: 


1. When the source charge distribution does not extend to ©, the 
potential can be defined as the work per unit charge done to bring a 


test charge, q,, from infinite separation to the point, r. 


The work done by a force F between r, and rz is definedas f Fd, 


so the potential is " 
> 
r oo 
s > +> > > 
V(r) = q E-d& = -q J E> df (3) 
oo > 
r 


Even when the source charge distribution does extend to infinity, 
some convenient position, ro, can be chosen, and potentials can be 
measured from that point. 


> 


— 


> > 
Vir) = q FE = de (4) 


> 
Fo 


2. The potential can also be introduced from its other relationship to 
the field, E, namely 


> 


E(t) = -grad (V(r) (5) 


The gradient, called grad, is just a vector derivative. 


These two definitions of potential are essentially equivalent. We will not 
do much with the second definition for the potential, V(r), but it is the 
more general of the two forms. This second definition is also the one used 
Most in upper levels of e & m study. One very useful result of this second 
definition is that the electrostatic field, F(r), at the point, r, points in the 
direction of steepest decrease of the electrostatic potential, Vir), at that 
point r. The force on a test charge at r always point directly down the 
potential hill. 





The first definition points out most clearly that, unlike the E field, which 
is a vector function of position, the electrostatic potential, Vir), is a scalar 
function of position. E(r) has both magnitude and direction at each point, 
r, whereas the potential, V(r), has only a value at each point, r. The fact 
that V is a scalar makes it very easy to use, which is one of the reasons the 
potential appears so often. 


Again, textbooks derive the potentials for various simple charge distribu- 
tions. A table follows: 


Charge Distribution 


Point Charge, q 


Potential 


ky a/r 






A = the line charge density 


Line Charge, A -2k, Aln(r/r,) 


Practice Exercise: Show that the potential due to a point charge follows 
from the integral definition of potential and the earlier equation for the 
electrostatic field due to a point charge. 


Analytical methods which solve electrostatics problems are generally 
limited to those which can be integrated easily. Numerical, computer- 
based solutions are not so limited. For a general charge distribution, you 
can break the distribution up into pieces (which can be treated as point 
charges or tine charges or whatever) and then just add up the fields or 
potentials due to each piece. The computer-based methods are very general 
and yet are based only on the simple equations for fields and potentials 
(those quoted above}. 


VISUALIZING ELECTROSTATIC FIELDS 


There are several ways to visualize electrostatic fields, E The first way uses 
the electrostatic potential, V. The starting point is a plot of the potential 
everywhere in space. (Everywhere means at lots of points throughout some 
region of interest.) You then connect all the points having the same value 
of potential. The resulting surfaces in three dimensions (or contours in two 
dimensions) are called equipotentials. No work is done on or by a test 
charge moved around on any equipotential surface. The forces on any test 
charge are always perpendicular to equipotential surfaces. Because the 
electrostatic field is the force per unit charge on a test charge, Q2,the E 
field is always perpendicular to equipotential surfaces. Thus you can 
visualize what happens to a test charge (in some field set up by given 
source charges) by finding the equipotentials in the field. You can also 
(crudely) find the electrostatic field, E, at any point, r, by discovering in 
which direction away from T the potential decreases fastest. The greater 
the rate of change of potential, the stronger is the electric field, E. 


All you need to be able to do to use this method is to be able to calculate 
the potential everywhere due to any given source charge distribution. 
However, the equation for an equipotential, is difficult to derive even for 
fields as simple as that due to three point charges. With the computer you 
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just add up the potentials due to whatever set of source charges you have: 
for three point charges, q,, 42, q3 at ri, r5 and r3 and for any point, r, 
you sum Kg (q,/lr -F | + q,/|? -r,| + q3/|r -73 |). For N charges, you add 
up N such terms. If the source distribution is some mixture of point 
charges, line charges and whatever else, you must remember to put the 
right form of the potential in each term. The method is easy and very 
general. A block diagram of this strategy for potentials is shown below: 







INPUT 


Charges and Their Positions 





> 


Set Range of r 


. For Each r, Calculate V (r) 


A second method is occasionally used to visualize electrostatic fields. In 
this method, the vectors E(r) are calculated at lots of points over the 
whole region of interest. Arrows representing the various fields, E(r), at 
the various points will be short where IE] is smal! and long where \E| is 
large. You can then visualize the E field in terms of constant IE] surfaces 
(or contours in two dimensions). Although this is a valid way to visualize 
fields, it is not very common. This method is illustrated by a problem at 
the end of this chapter. 









The third and most common way to visualize electrostatic field patterns 
uses field lines. One of the main points of this chapter is to allow you to 
use simple computer programs to generate field line patterns for numbers 
of different source charge distributions. 


The two-part definition of a field line is deceptively simple: 


—_ 
1. At every point, r, the electrostatic field, E, is tangent to (that is, is 
locally parallel to) the field line. 


2. The number of field lines passing through a (unit, perpendicular) 
> 
cross sectional area at r is proportional to |E({r)|. 
This second part of the definition simply means that the field is strong 
where field lines bunch together and weak where lines are spread apart. 


Notice a few points about field lines: 


=> 
a. {E| is usually not constant along a field line. Any line usually goes 
—> 
through regions where E is weak and regions where E is strong. 





b. Field lines start on positive charges and end on negative charges. If 
the total net charge of a charge distribution is not zero (that is, if 
there is more positive charge than negative, or vice versa) then some 
field lines may begin or end at infinity, since the lines will eNep on 
searching for some charge on which to terminate. 


c. The number of fines that start on a +2 charge is twice the number 
that start on a +1 charge. Similarly, twice as many lines will end ona 
~2 charge as end on a -1 charge. You can choose arbitrarily how 
many lines start on a +1 charge (note that the second part of the 
field tine definition says “‘proportional to’’), but having chosen that 
number, all else follows. 


At first sight the definition of field line seems very simple: after all, there 
are only two defining characteristics. At second glance, the definition 
seems to make field lines very complicated. Analytical solutions for field 
lines can only be performed for very simple source charge distributions, 
and even then the expressions look very complicated. Using a computer 
the field lines can be mapped using just the definition and a very short 
program, 


Consider the following way to trace out a field line. Suppose you know 
some point, r= (x,y,z), on a field line. (Any point lies on some field line, 
so you define the field line you are going to trace by choosing the first 
point.) Suppose further that you want to take a step of size As along that 
field line. You calculate the field at r, E(r), and then use the first defining 
property (the line is parallel to E at r) to calculate 


—-> —> 

Ax = As E,/|E| 
> 

Ay = As Ey/|E| 


~> 
Az = As E,/|E| (6) 


: — —_> > —-> - =a oe >. 
These equations say that As = Ax + Ay + Az is parallel to E = Ey Ey tle 
which is equivalent to the first defining property. 


Practice Exercise: When e =0, show that the triangle Ax, Ay, As is similar 
to the triangle E,, E y? E 


The next point along the field line is then (x+Ax, ytAy, z+Az), and you 
can repeat the process at this new point. You walk along the field line 
step-by-step. Such a method of solution is called “algorithmic” (because 
the solution is an algorithm or procedure) or “‘iterative’’ (because you 
iterate or repeat the same steps over and over). Notice that the procedure 
breaks down if IE | = 0 at any point on the field line. Such methods are 
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easy on computers: A block diagram of this procedure is shown below: 


INPUT 
Charges and their positions; 
Step size, As 


. > 
Choose point, r, on 


field line 


Galculate the total field, 


E,atr 


Calculate Ax = As E,/E 
Ay = As EV/E 
Az = As E>/E 


Calculate new r = 
(x + Ax, y + Ay, z + Az) 





Notice that you could perform the calculation yourself if you had the 
patience; there is nothing magic about the method at all. 


What about the second part of. the field line definition? How do you 
ensure that your pattern will obey that property, too? The answer lies in 
how you choose to start field lines. Close enough to any charge, the field 
lines wil! be radial since all the other charges are so far away their effects 
are negligible. Close to each positive charge you can choose the number of 
lines (q times the number you choose to start on a +1 charge) to start 
radially and at equal angles. Once the lines are started correctly, the whole 
pattern will remain correct. More Jines will bend around so that they go . 
through regions of strong field, and fewer lines will bend around so as to 
go through weak field regions. This simple method produces results 
automatically. 


The same sort of procedure can be used to trace out equipotentials. Since 
equipotentials are always perpendicular to field lines, we move perpen- 
dicularly to the field, E, at each point (rather than parallel to E). in two 
dimensions (or symmetrical three-dimensional systems such as those we 
consider below) this procedure is particularly easy. The steps Ax and Ay, 
parallel to x and y, are given by Ax = -As EY /|E| and Ay = As E,/|E}. 





FIELD LINES AND EQUIPOTENTIALS WITH’ THE 
COMPUTER 


Let us do an example of this procedure which is similar to the problems at 
the end of the chapter. We will choose a relatively simple example: two 
opposite charges (+1 and-1) at (+.5, 0) and (-.5, 0) respectively. Far away 
from the charges (\7[>>1) this pattern reduces to a dipole pattern. (This is 
actually one of the few problems which can be solved analytically. Hence 
this example can also serve to verify the numerical method; this verifica- 
tion is discussed in a problem.) 


Since the pattern is symmetric under rotations around the x-axis, we can 
limit our discussion to the x-y plane. That is to say, the z direction is just 
like the y direction, so the full three-dimensional picture is just the one we 
will produce rotated around the x-axis. Two-dimensional pictures are 
easier to put on paper; the method, of course, works just as well in three 
dimensions. Most of the physical intuition can be derived from two- 
dimensional pictures. 


Example 1. Find the potential for the above charge distribution. Consider 
-5<x<5 and -5<y<5. 


The potential is just V = ko (+1/[distance from (+.5,0)]) + kg (-1/[distance 
from (-.5, 0)}). One program which calculates and prints-out the values of 
the potential is listed below: 


POTENT 


1@0 PRINT *€XsyY) OF LEFT-BOTTOM & eet Specify Region 
11@ INPUT X8,Y8,X9,Y9 of interest 


122 
132 
148 
150 
168 
17@ 
188 
198 
208 
218 
228 
238 
242 
258 
268 
27 @ 
288 
29@ 
328 


Practice Exercise: Modify the program so that it computes the potential 
for two equal, positive charges at (+.5, 0), {-.5, QO). 


LET K@= [-— Set units 

LET Q=+ ! — Magnitude of changes 

LET N@29 — #of points across x and y 

LET X73(X9<X8) /N@ — Ax 

LET Y7=CY9-Y8) /N@ — Ay 

FOR J=N@ TO @ STEP <1 — Step downy 

LET YO=Y8+JI*Y7 — Each y in turn 

FOR 128 TO N@ — Step across x 

LET X@2X&+1*X7— Each x in turn 

LET R=SQRC(X@~--5)*(XO-- 5)4+YO*xYOQ) — Distance from (+.5,0) 
LET Ri®SQRCCXG+-5)*CXO+2 5) +Y OXY) — Distance from {-.5,0} 
LET V=K8*Q/R-K@*Q/R1 — Potential 

PRINT Vz — Print potential 


NEXT I 
PRINT 
NEXT J 
PRINT 


3O0TO 1@@ — Return for new region 


END 
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A RUN of this program produces the following output. Notice that each 
value of y produces two lines of printing; the second line gives the values 
another page width to the right of the first line. 





RUN 
POTENT 


(XeY) OF LEFT-BOTTOM & RIGHT-TOP? 


20 5 0 50 2 55 eS 


-~t- 18557 -.971852 -+ 746188 -- 46418! ~+« 156889 
« 156889 4641 «746108 °971852 1-!8557 
-1le 63942 ~1244181 -1-08265 -+ 656786 --218 68 
«21868 « 656706 1688265 12-4418 1 1.63942 
72+ 63648 -2+ 2687 3 -1+68832 --928048 ~ +298 828 
« 298028 «9288648 1.68832 2+ 26873 2+ 63648 
-5.81361 -3-88657 -2-34283 - 1.22887 ~-e382654 
© 382654 1-228 87 22342839 3288657 5°81361 
“17-8615 ~6-92764 -3-08319 ~1-46436 - «441558 
« 441558 1+ 46436 3- 88319 692763 17-8815 
-17-8815 -6+927B4 -3+@8319 -1+46436 - -441558 
«441558 1+ 46436 3¢ 88319 6292783 17-8815 
-5-81361 -3-88657 2234283 -1-22887 --382653 
° 382683 1+ 228 87 234283 3°88657 $-8@1361 
-2+ 63648 ~2- 26873 ~1. 68832 ~»9 280648 - + 298828 
© 298828 920048 le 60832 2+26873 2+ 63 648 
-1-63942 ~1-44181 * -1.88265 ~- 656766 - +218 68 
- 21868 » 656786 1- 88265 12644181 1-639 42 
~1le18557 -~097 1852 -»- 746188 -2464161 -»+ 156889 
» 156889 4641 «746128 971852 1-1@6557 
(XsY) OF LEFT-BOTTOM & RIGHT-TOP? 
2058s 141 
8 787 228E-82 e 152839 ©218173 © 271447 
© 318716 © 335683 © 34721 ©3477 2339727 
a «18363 « 286662 285202 °351774 
- 398157 © 423933 «431343 “424242 486995 
@ © 138997 + 268 648 ©3799 88 » 465189 
«5197 62 «543989 «542368 521774 ° 489 679 
a ¢ 1899 37 ° 367195 518171 «629794 
+ 69 38 29 711685 » 69 87 67 © 646572 °598792 
@ © 263886 °512489 ° 725627 °87843 
©9527 23 958207 892262 *-8@6114 712763 
@ - 369911 728777 1+-8@4792 1.27458 
163595 le 38575 1.16862 1288755 85562 
6 515818 1.84749 21-5691) 1-9 6872 
285579 1.859 12 165474 1.25226 1683331 
QB ° 69478 1 1.48774 2244652 323349 6 
3+ 43859 2-758 2+ 84289 1«52167 1-l6815 
@ «862522 1.974 3-8 8283 6299827 
7.187 68 4-¢13982 2.56285 1.75478 1.28753 
@ 935665 2.21538 48 16-9412 
17-8526 5° 14286 2.81739 1-85143 1+33333 


(X,Y) OF LEFT-BOTTOM & RIGHT-TOP? 
? 


The largest positive numbers for the potential! lie near the positive charge, 
the largest negative values of potential lie near the negative charge. 


Sometimes plotting the potential helps you to see the pattern. You can 
plot characters on the terminal which represent the size of the potential at 
each {x,y} point on the plot. One way to produce such a plot involves 
several changes (shown below) in the potential program, POTENT, listed 


on page 8. The characters represent the size of the potential; the position 
of each character is the point (x,y). Terminal plotting is sometimes crude 
since on most terminals the characters can only appear at definite places 


on the paper. 


GET~-POTENT 


10 DIM A$(36) 

28 LET A$#"@123456789 ABCDEFGHI JKLMNOPQRSTUVWXYZ" 
[12 PRINT "MINe & MAXe V7?" 

114 INPUT V8, VW 

14@ LET N@=36 

168 LET Y7=16*(Y9-Y8) /( 6*#NB) 

178 FOR J#INTC. 6*#NQ+-5) TO @ STEP -1 
249 LET I9S1+INTC 35% V- VB) /(V9-V8B)+65) 
241 IF I9>"81 THEN 243 

242 LET [981 

243 IF 19<836 THEN 245 

244 LET 19236 

245 PRINT ASC(I9,19)3 


A RUN of this modified program looks like: 


RUN 
POTENT 


(XsY) OF LEFT-BOTTOM & RIGHT-TOP? 
Im 0 854-0 852-2 45,205 

MIN & MAX V7? 

27-1800 
UUUUUUTTTTTTSSSSSSSTTTTTTUUUUUU 
UUUTTTTTSSSSSSRRRSSSSSSTTTTUUUU 
UUTTTTSSSRRRRQQQQQRRRRSSSTTTTUU 
UTTTSSSRRAQGPPPPPPPAQARRSSSTTTU 
TTTSSRRQQPPONNNNNNN OPPQGRRSSTTT 
TTSSRRQPOONMLKKKKKLMNOOPQRRSSTT 
TSSRRQPONMKJHGFEFGHJLMN OPQRRSST 
TSSRQPONMKI FC8 S458 CFIKMNOPQRSST 
SSRRQPOMLIFBS@@eG8SBFILNOPQARRSS 
SSRROPOMKIES@GA2ZAG29 EIKMOPQRRSS 
SSRRQPOMLI FBES@@egagSBFILNOPQRRSS 
TSSRQOPONMKIFC85458 CFIKMN OPQRSST 
TSSRRQPONMKJHGFEFGHJLMN OPQARRSST 
TTSSRRQPOONMLKKKKKLMN OOPQRRSSTT 
TTTSSRRQQ@PPONNNNNNN OPPQORRSSTTT 
UTTTSSSRRQQQPPPPPPPAQOARRSSSTTTU 
UUTTTTSSSRRRRQQQQAQRRRASSSTTTTUV 
UUUTTTTTSSSSSSRRRSSSSSSTTTTUVUY 
UYUUUUUUTTTTTTSSSSSSSTTTTTTUUUUUU 


(X,Y) OF LEFT-BOTTOM & RIGHT-TOP? 
Fo Sa 0 $92 See S 

MIN & MAX V? 

2-252 25 
HHHHHHHHHHHHHHHITIILTIIIIIIIIII1 
HHHHHHHHHHHHHHHIIILTIIIIIIIIIIII! 
GGGGHHHHHHHHHHHIIIITITIIILLI ddd 
GGGGGGHHHHHHHHHITITITIIIIII I JJ 
GGGGGGGGHHHHHHHIIITITIIJI JI JIU 
FFF FGGGGGHHHHHHIIILIITIJVJJJJKKKK 
EEEFFFGGGGHHHHHIITITIJVJJJKKKLLL 
CCDEEFFGGGHHHHHIIIITIJJJKKLLMNN 
67 ACDEFFGGGHHHHITIIIJJJKKLMNPST 
208 BDEF FGGGHHHHIIIIIJJJKKLMORZZ 
67 ACDEFFGGGHHHHIIIITIJJJKKLMNPST 
CCDEEFFGGGHHHHHIITIIIIJJJKKLLMNN 
EEEFFFGGGGHHHHHIIIITIJJJJKKKLLL 
FFF FGGGGGHHHHHHITIITILIJJSJJKKKK 
GGGGGGGGHHHHHHHILTIIITITIJJJJJJJJ 
GGGGGGHHHHHHHHHIIIIIIIITIiddJJJJ 
GGGGHHHHHHHHHHHITITIIITIITILIIS dd 
HHHHHHHHHHHHHHHITIIITIIILIIIIIII 
HHHHHHHHHHHHHHHITIIIILZIIIIIIIII 


(XsY) OF LEFT-BOTTOM & RIGHT-TOP? 
? 


(The character @ represents V< = the minimum potential; the character Z 


represents V> = the maximum potential.) 


Example 2. Sketch three equipotentials (approximately) on a plot of 


potentials. 
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Using the terminal plots you can connect the points having approximately 
equal values of potential. Points having the same character are connected 
(but not @ or Z); sometimes it is necessary to interpolate between char- 
acters (for example to find where 6 would come between printed 5 and 7 
characters). The breadth of a region having a given character is broadest 
when the potential is changing slowly. 


Practice Exercise: On the second terminal plot, find the boundary be- 
tween the | and J equipotentials. What is the shape of this equipotential off 
the plot? 


Practice Exercise: On the first terminal plot, find the K equipotential. 


Practice Exercise: Find the contour on which the value of the potential is 
actually zero on the first terminal plot. Explain its behavior. (/t is possible 
to have the computer hunt through and interpolate between points to find 
points ona given equipotential.} 


Example 3. Sketch three field lines (approximately) on a plot of 
potentials. 


You can find field lines by moving in the direction of the fastest decrease 
in potential at each point. Using terminal plots the field lines are crude, 
but you can still get a feeling for this behavior. 


Practice Exercise: Starting near the lower left corner of the first terminal 
plot, follow a field line inward. 


Practice Exercise: On the second terminal plot, follow a field line which 
leaves the positive charge at 135° to the x-axis. 


Example 4. Using an algorithmic method, compute three field lines for 
this charge distribution. 


You can use a program which implements the algorithmic strategy 
discussed above. One program which does this is the following: 





Initialization 







Main Intergation Loop 


All calculations of the distance from each charge and of the fields use an 
estimate of the position one-half step ahead of the present point on the 
line. This half-step method is quite accurate even with a relatively large 
step size {and hence relatively few iterations to complete a line). The 
question of accuracy and the idea of a half-step method are discussed 


LET NG@=2 — #of charges 

FOR J=1 TO N@ 

READ X€J1,YCJ1,0€J] Positions and 
NEXT J Values of Charges 
DATA «520s lsa-25s 8,71 

LET K@29.E+@9 — Units 

LET D=.@5 — Step size, As 

PRINT "*€XsY) STARTING PT- ON FIELD LINE?" 
INPUT X5,Y5 

LET X@=XS5 — x along line 

LET Y@=Y5 -—- yatong line 

LET X1=@-— Ax 

LET Y1i=@-Ay 

LET S@=.-5 — Print out every S@ along line 

LET S=@-— Distance S along line 

PRINT Xs ¥ * 

LET E1204 

LET E22@ 

FOR J=1 TO N@ —- Add up fields for each charge 
LET R3=(XG+X1/2-KCI7)*CXB4+X1/2-XC57) |e 
LET R3=R34+ (CY G+Y1/2-YLd 2) CY RtY 1 72-¥CU7)) 
LET R3=R3t t.5—};)3 

LET E=K@*Q(J)/R3—lEIAr 

LET EI=E1+(X8+X1/2-XCJ1])*E- E, 
LET EQ=E2+ (YOY 1/2-YCIII*E- EY 
NEXT J = 

LET E@=SQRCEI*E1]+E2*E2)~— iF 

LET X1=D*E1/E@-— Ax 

LET X@=X@+X 1 — New x 

LET Y1=D*E2/EO—- Ay 

LET Y@=YO+*yY I- Newy 


LET S=S+D 
LET S=@ Lines 


PRINT X@,Y2 

FOR J=1 TO N@ 

IF ABS(CX@-XCJ1I)+ABS(YB-YC(J3)<-9*"D THEN 492 
NEXT J 

IF ABSC(X@)+ABS(Y@)<1@ THEN 268 - Far off page 
PRINT X@sYQ@ — Print out last point 

PRINT 

GOTO 172 — Aeturn for new line 

END 


further in Appendix B. 


A RUN of this program looks like the following: 


EMAP 


(XsY) STARTING PT. ON FIELD LINE? 


2e555-85 
Xx 
«8797 27 
le 13562 
le 2977 
1635226 
le 28674 
1+ @9879 
«7 62383 
321695 
-- 173584 
~+ 638 683 
-1-8@@482 
~1e24142 
-1.34563 
“1.32623 
~121954 
=+9 66379 
™e 655836 
-2515865 


Y 
425224 
853816 
1-32582 
1682174 
2- 31608 
2077418 
3+ 14798 
3637725 
32 489 38 
3223573 
2e89915 
20 46887 
1.97332 
1247485 
©99 326 1 
549795 
° 158764 
815919 


(XsY) STARTING PT- ON FIELD LINE? 


20555 2-8! 

x 
1¢683884 
1.52333 
2e 88277 
2e AT 694 
2094573 


Y 
© 114918 
© 238369 
«38O176 
«5387 25 
712554 
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3¢ 489 87 °9 80448 
3-8 6689 12-1814 

4.31915 1.31457 
4-7 6582 1653924 
5+ 28685 1277478 
5« 64222 2e 82864 
6-87 188 2- 27634 
6249579 2054145 
6.91394 2.81558 
7+ 87958 2e9 27 67 





(X.Y) STARTING PT. ON FIELD LINE? 
? 


The field tines start at the positive charge and end at the negative charge. 
Some lines are much longer than others; the longer ones bend around from 
the back side of the positive charge to the back side of the negative charge. 
The line leaving the back side of the positive charge directly away from the 
negative charge must go through infinity to get back to the negative 
charge. 


Practice Exercise: Modify the program to find field lines for four point 
charges: +1 at (+1, +1), -1 at (-1, +1), #1 at (-1,-1) and -1 at (+1, -7). 


Plotting can be helpful in this exercise. If you have a plotter system or a 
plotting terminal, use it; if not, terminal plotting can serve. The changes 
required to use terminal plotting with EMAP refer to subroutines in 
Appendix A and are as follows: 


GET-EMAP 

18 READ X8sX9sY8sY92Z8 
28 DATA <-5s52-5s 5s 1 

38 GOSUB 98988 

285 GOSUB 9168 

258 

448 GOSUB 9188 

498 GOSUB 91608 

ag92 GOSUB 9288 


APP-TTYPLO 


RUN 
EMAP 


(XsY) STARTING PT. ON FIELD LINE? 

72-55,-05 

MAX Y® 5 
XXXXXXXXXKXKXXXXXAXKXKX XAXXXKXX XX KXXXKX AAKXXXXKKXKAXKXXKKKX 


SSeS KKK KKK KR KR KH 


XXXXKKXXXKXXKXKKXKXAXXKKXXKXXKXXKXKXKXKAKKKAAKKAXA KK AXKXKKAKKAKKX 
MIN Y#-5 MIN X#-5 MAX X= 5 
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(XeVY) STARTING PT. ON FIELD LINE? 
Ve 55, ese a3 


MAX Y= 5 
XXXXXXXXXXXKXXKXXXXXXXXXXXKXXXXKXKXKKAXKXXKKXKKAKKXKARKAK 


} 1 
l 


Midi dee Keddie ee Kedah KK KKK KKM KKK KKK 


1 
XXXXXXXXXKXKXX XXX XXX KKK XX XXX KX XKXXXKXKXKX KX KKK XXX KXKXAXKKXAXXK 
MIN Y=-S MIN X=-5 MAX X= 5 


(XsY) STARTING PT. ON FIELD LINE? 


XXXXXXXXXXXXXXXXXXXXXAXXKXKXKKXXXKKAXAKKKK KKK KAKKKKXKAKA 


Khe d eee KKK KKK KKK KKK KKK 


XXXXXXXXXXKXXKKXKXXXKXXXKXKXKXKXXKXKAXAXK AK KK KKAKKXKXKKKAKKX 
MIN Ys-5 MIN X==<5 MAX X= 5 


(XsY) STARTING PTe ON FIELD LINE? 


= 


Sometimes it is useful to have several lines on one plot. The following 
changes in the original field line program, EMAP, listed on page 11, 
accomplish this. 
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GET-EMAP 

18 READ X8s5X9sY8sY95Z20 
29 DATA -5,55-525,08 

38 GOSUB 9882 

182 IF X52999 THEN 515 
183 LET Z@=Z8+1 

285 GOSUSB 912828 

2528 

442 GOSUB 9189 

49@ GOSUB 9188 

515 GOSUB 9286 


APP-TTY PLO 


A RUN of this modified program is shown below: 


RUN 

EMAP 

(XsY) STARTING PT. ON FIELD LINE? 
70552285 

(XeY) STARTING PTe ON FIELD LINE? 
2705527283 

(Xs¥) STARTING PTe ON FIELD LINE? 

204526 

(XsY) STARTING PT- ON FIELD LINE? 

799958 

MAX Y= 5 
XXXXXXXXXXKXKXXXXXXXXAKKXAXKXAXKKXXAXXKXXXXKXAXAXKKXXXXKXXAXK 
Y 

Y 

Y 

Y 

Y 

Y 1 i 4 

Y I l 

Y 1 

© i 1 

Y i 

Y 1 

Y l i 

Y t 1 

Y 1 

Y 1 

mt 23 3 32 

Y 2 2 

Y 2 2 

Y 2 2 

Y 2 

* 2 2 

Y 2 

Ys 2 2 
Y 2 
¥ 2 

Y 2 2 
Y 2 
Y 2 

Y 2 2 
Y 2 
Bg 2 
XXXXXXXXKKXXXXK XX XXX KKXXX XXX XXXXXXKXKXXAXKAKXKAXAAXAXAXAAXX 
MIN Y#e-S MIN X#=#<-5 MAX X= 5 


(The terminal plotting programming prints the last character plotted at 
each position.) 


Finally, it is possible to modify the field line program to follow equi- 
potential contours instead. All you do is move perpendicularly to the field 
E at each point (instead of parallel). The modified program is listed below, 
followed by a run and a terminal pilot. 








ELECTROSTATICS O 15 


VM AP 
1@0 LET N@=2 
11@ FOR J21] TO N@ 
120 READ X(J1],YC€J),QCJ) 
138 NEXT J 
1498 DATA ©52Gs les 2°54 8s-1 
150 LET KO0=9-E+99 
168 LET D=.85 
176 PRINT "(Xs Y) STARTING PT» ON FIELD LINE?" 
186 INPUT X5.Y5 
1986 LET X@=X5 
2080 LET Y@=Y5 
218 LET X1=@ 
226 LET Y138 
2380 LET S@=.5 
248 LET S=8 
25@ PRINT ‘°X's ‘ty 
260 LET El=@ 
27@ LET Ee=e 
286 FOR J=i TO N@ 
296 LET R3=(XO4+X1/2-XCSI]>*¢(X8+X 1 /2-XCSI)) 
308 LET R3=R3+CY84Y1/2-YCSI)*«(YO+Y1/2-YC5)) 
310 LET R3FRIt1-5 
320 LET E=K@*QLJ3/R3 
33@ LET EL=EL+(X8+X1/2-XCUII#*E 
340 LET Een E24+(YO+Y1/2-YCJI*E 
358 NEXT J 
368 LET E@2SQRCEI*E]+E2*eE2) 
378 LET X12-D*E2/E@ 
389 LET X82XO+X! 
39@ LET Y1=D*E1/EG 
4@6 LET Y8=YO6+Y1 
416 LET S=S+D 
426 IF S<S@ THEN 4580 
430 LET S36 
440 PRINT X8sY@ 
450 IF ABS(X@-X5)+ABS(YO-Y5)<+9*D THEN 478 
460 IF ABS(CX@)+ABS(Y@)<1!6 THEN 268 
47@ PRINT XOsY@ 
48@ PRINT 
498 GOTO 178 
588 END 
VM AP 
(XesY) STARTING PT. ON FIELD LINE? 
263,98 
x x 
« 655171 -- 166054 
© $444) 219985 
° 293841 3°87 172E-82 
(XsY) STARTING PT. ON FIELD LINE? 
72158 
x Y 
° 263284 -- 459612 
719458 =~» 618559 
1612219 -+ 339 661 
1-28455 © 148735 
921518 «536615 
« 449887 «58 67 65 
« 127882 ©214384 


9-85179E-982 


(XeY) STARTING 
7-e158 

x 

“e 263284 

-«7 19458 
-1-12219 
-1.-28455S 
~0921519 

“~e 440887 

-- 127882 

~9 -BSITSE-82 


(XseY) STARTING 
? 


1-68 651E-82 


PT» ON FIELD LINE? 
Y 
-- 459 612 
~» 618559 
-« 339 661 
e 148735 
536616 
«589765 
«214384 
1.68 656E-02 


PT. ON FIELD LINE? 
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GET=-VMAP 

1@ READ X8,X9,Y8s,Y9.Z208 
2@ DATA -3, 34735. 358 

38 GOSUB 9866 

18S IF X5=999 THEN 495 
186 LET Z@=Z@+1 

205 GOSUB 9188 

258 

448 GOSUB 91828 

47@ GOSUB 91828 

495 GOSUB 923¢ 


APP-TTYPLO 
RUN 
VMAP 
(XY) STARTING PT- ON FIELD LINE? 
723,80 
(X.Y) STARTING PT. ON FIELD LINE? 
201,90 
(XY) STARTING PT- ON FIELD LINE? 
7-2 1,2 
(XY) STARTING PT. ON FIELD LINE? 
799958 
MAX Y= 3 
XXXKXXXKKXXXXXXKAXKK KKXXXKAKKKXKKXKKA KK KX KK KXAKXK AK XK K KAR 
7; 
Y. 
v 
x 
¥ 
&¢ 
Y 
Y 
Y 
Y 
Y 
om: 
Y 3 3 2 2 
r 
~ 3 3 2 1 2 
Y 3 211 
* 1 
Y 3 3 2 2 
Y 3 2 
Y 
Y 
Y 
¥. 
¥ 
Y 
Y 
Y 
Y 
¥ 
+ 
Y 


XXXKXKXXXXAAXXKARKXXXKXXXAXXKKKXXKXKXAKAXKKAKAKKX KX KKKXKXKKKKK 
MIN Y#-3 MIN X==-3 MAX X23 3 


Practice Exercise: Annotate the program with comments by each line or 
group of lines explaining what calculation the line or group perform. You 
may use the annotations on the earlier program as a guide. 


Practice Exercise: Modify either of the algorithmic programs so that the 
program follows either a field line or an equipotential. Allow the user to 
decide which type of contour to follow each time he starts a new contour. 


CONCLUSION 


The methods discussed here are very general. Field lines and equipotentials 
can be found for any source charge distribution. By writing and running 
your own programs for several different charge distributions, you will gain 
a very good feeling for the behavior of field lines and equipotentials. Field 
line mapping is fundamental to every part of electricity and magnetism, so 
it will be helpful to gain this intuition now. 





CHAPTER ONE EXERCISES 


1. Consider the line charge distribution: +3 at (+.5,0) and -1 at (-.5,0). 
(This potential is symmetric around the x-axis so we only need to 
view it in the x-y plane). 


a) 


b) 


c) 


d) 


e) 


Find the potential for this charge distribution. Use the char- 
acter plotting method to observe the potential. (You will need 
to print the values in the region first in order to know the 
largest and smallest values of the potential to plot.) Consider 
the regions: -.55 <x <-.45, -.05 <y <.05,; +.45 <x <.55, 
-.05 <y <.05; and -5 <x <5,-5<y<&5. 


Sketch several equipotentials on each of the character plots. 
Since the change in potential is linear with the character set, if 
you sketch the equipotentials for every third character (eg., 2, 
5, 8 etc.), then the contour lines are close together when the 
change in the potential is large (and the field is strong). 


Sketch four field lines on the character plots. 


Using an algorithmic method, compute four field lines for the 
large region considered in (a). 


Using an algorithmic method, compute four equipotential lines 
for the large region considered in (a). 


2. Consider the line charge distribution A = +1 at (+.5, -.25), X= +1 at 


(-.5, -.25) and A = +1 at (0, +.62), an equilateral triangle. The line 
charges are perpendicular to the x-y plane. 


a) 


b) 


c} 


d) 


e) 


f) 


Find the potential for this charge distribution. Use the char- 
acter plotting method. Consider the region -2<x <2,-2<y 
<2. 


Sketch several equipotentials on the character plot. 


Sketch three field lines starting at each charge on the character 
plot. 


Using an algorithmic method, compute four field lines starting 
at each charge of this charge distribution. 


Using a algorithmic method, compute several equipotential 
lines. Pick equipotentials close to each charge and then far 
away from all of them. 


How would you compare the potential near each charge and 
far away from all of them to a single line charge distribution? 
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3. Consider the line charge distribution of \ = +1 at (.5, .5);A =-1 at 
(.5, -5); A = +1 at (-.5, -.5); and A = -1 at (-.5, +.5). The line 
charges are perpendicular to the x-y plane. 


a) Find the potential for this charge distribution. Use the char- 
acter plotting method. Consider the region -2 <x <2, -2< 
y<2. 


b) Sketch several equipotentials on the character plot. 


c} Sketch three field lines starting at each charge on the character 
plot. 


d) Using an algorithmic method, compute four field lines for this 
charge distribution. 


e) Using an algorithmic method, compute several equipotential 
lines. Pick equipotentials close to each charge and then far 
away from all of them, 


F) Locate the planes on which the potential is zero. 


4. Modify Exercise 3 so that each of the line charges has A = +1. Sketch 
your idea of the equipotential and the field lines before you calcu- 
late them, 


a} Using an algorithmic method, compute four field lines starting 
at each positive charge for this charge distribution. 


b) Using an algorithmic method, compute equipotential lines for 
the distribution. 


c) How would you compare the potential near each charge and 
far away from all of them to a single charge distribution. 


5. Let us examine the idea that a set of point charges can be used to 
approximate any field distribution. To do this, let us approximate a 
line charge by a set of point charges. Consider the point charge 
distribution with nine +1 charges from (-2,0) to (2,0) at .5 intervals. 


a) Compute the field line distribution for this charge distribution. 
b) Compute the equipotential lines for this distribution. 


c) Since you know the analytical expression for the field and 
potential from a line charge, compare your point charge ap- 
proximation to it. Where does the point charge approximation 
to field lines and potential differ from the line charge? (There 
are three regions where the two differ most). 


d) How would you make a better approximation to the line 
charge? 


Consider the line charge distribution of \ = +1 at (0,0); A = -1 at 
(-1,0); A = -1 at (+1,0); A = +1 at (-2,0); A= +1 at (2,0); A =-1 at 
(-3,0); X= -1 at (3,0); A = +1 at (-4,0); and A = +1 at (4,0). 


a) Compute the field line distribution for this charge distribution 
of -.5< x <.5,-1<y<7. 


b) Compute the equipotential lines for this distribution. 
Cc) On what surfaces is the potential zero? 


In Exercises 3 and 6, and in the dipole example, the charge distribu- 
tion leads to planes of zero potential. These problems illustrate the 
image phenomena. The field pattern of these charge distributions is 
the same as that of conducting sheets at zero potential and the 
charge distribution on one side. For example, Exercise 3 can be 
thought of as +1 at (.5, .5);-1 at (.5, -.5) and a conducting plane at 
zero potential along the y-axis. The way to solve the conducting 
sheet problem is to set it up as if each charge has an image of 
opposite sign an equal distance on the other side of the plane. 
Thus, Exercise 3 is the solution to the two charges and the con- 
ducting plane. 


a) Exercise 3 is also the solution to another charge distribution 
and a plane. What is it? 


b) — Exercise 3 can be seen as the solution to the charge distribu- 
tion for two conducting planes and a point charge. Where are 
the planes? 


c) How can Exercise 6 be restated in terms of conducting planes? 

d) Consider the line charge distribution \ = +2 at (.5,0), \=-17 at 
5, .5); X\=-17 at (.5, -.5) and a conducting plane along the yz 
axis. How would you solve this problem using the method of 
images? 

Consider the charge distribution caused by an infinite solid cylinder 

of charge centered on the z-axis with a radius of .5. This can be 

approximated by a set of line charges parallel to the z-axis. 

a) Compute the field line distribution for this charge distribution. 

b) Compute the equipotential lines for this distribution. 

Cc) What is the potential inside the cylinder? 

d) What should the field be according to Gauss’s law inside the 


cylinder? How does the potential and field outside of the 
cylinder compare to those of a single line charge? 


We have visualized fields by either equipotentials or field lines. We 
can also visualize fields by surfaces of constant field magnitude. 
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10. 


Consider the dipole charge distribution we used in the text. 


a) Print the magnitude of the field over the region used in the 
text. 


b) Convert the character printing program so that the characters 
represent magnitudes of fields. 


c} Draw in the equifield regions on the character plots. 
d) How do the equifield regions differ from equipotential lines? 


e) Is there any relationship between the constant field contours 
and the directions in which the fields point? 


The analytical solution of a dipole is: 
the potential V = pcos 6/r? 


and E, = 3psin @ cos 6/r° 
E,= p (3cos? 6 -1)/r° 


are the field components when the dipole p is centered at the origin 
and is parallel to the z direction. The angle 8 is measured with respect 
to the z direction, and r is measured from the center of the dipole. 


If the dipole is caused by two single charges on the z-axis, then [pj = 
qs where s is the distance between the charges (the positive charge is 
assumed to be in the positive z direction). 


a) Compare the program results obtained for the potential calcu- 
lated from two point charges with the analytical dipole ap- 
proximation. (Compare the values on the x- and z-axes at the 
distances 1, 5, 20.) 


b) Compare the calculated results for the field with the analytical 
approximation. To do this, see if the two solutions are parallel 
to each other along the x-axis at 1, 5, 20 and along the 45° 
line between the x- and z-axes at the distances 1, 5, 20 from 
the origin. 





CHAPTER TWO: MAGNETOSTATIC 
FIELDS 


INTRODUCTION 


Magnetic fields Bir) result from currents. (We will use the term magnetic 
field even though magnetic induction is strictly correct; the distinction is 
really only important when you deal with magnetic materials.) If the 
currents are steady (that is, do not change with time), the magnetic field 
is magnetostatic. The purpose of this chapter is to find how to map 
magnetostatic fields by means of flux lines. 


MAGNETIC FIELDS FOR VARIOUS CURRENT 
DISTRIBUTIONS 


Physics textbooks derive equations for the magnetic field due to simple 
current distribution, such as the long straight wire. The texts also intro- 
duce the Biot-Savart law (sometimes called the Biot !aw) which gives the 
magnetic field everywhere in space due to a short segment of current, 1d@. 
These are the only results for magnetic field that we will need in this 
chapter. The Biot-Savart law is sufficiently general and, when you use a 
computer, sufficiently easy to use that any current distribution can be 
handled. 


For a long straight wire the magnetic field is given by B = 2k, l/r 
(tangential around the wire with sense by the right-hand rule), where r is 
the perpendicular distance to the wire, and kg is a constant which 
determines the units. ky = Uo/4a = 10°” in mks units in which B is 
measured in Tesla; Ky = 1/c = 1/3x10'° in Gaussian (cgs) units in which 
B is measured in Gauss. In mixed units where currents are measured in 
amperes and magnetic field in Gauss, kg = .1. (One Tesla is 10* Gauss.) 


—> 
The Biot-Savart law states that the magnetic field, dB, due to a small 
element of current, 1d0, IS 


> > >>, 
dB = ko ! d& x r/Ir| (7) 


where r is the distance from ag to the point, (x, y, z), that the field, aB, is 
measured, This relationship allows you to compute the magnetic field due 
to any circuit by breaking the Circuit up into little pieces and then adding 
up all the dB‘s. This relationship for dB is an inverse square law but has the 
complication of the vector cross product in the numerator. The cross 
product (in Cartesian coordinates especially) is quite easy to calculate ona 
computer. 
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MAPPING FIELD LINES 


There are two ways to visualize magnetic field patterns. The first is not 
often used but is easy to understand. In this first method, the magnetic 
field, B, is computed at a large number of points throughout the region of 
interest. The entire field is then visualized in terms of contours (or surfaces 
in three dimensions) of constant BI. This method is not in common use 
but is considered further in a problem at the end of the chapter. 


The second (and very common) way to visualize magnetic field patterns 
uses flux lines. The definition of a flux line has two parts: (1) the 
magnetic field, Bir), is tangential to the flux line at each point, r, 
and (2) the number of flux lines passing through a (unit, perpendicular) 
cross section placed at any point, r, is proportional to the magnitude of 
the magnetic field, |B(r)|, at that point. 


To calculate the shapes of flux lines analytically (the result then being an 
equation for the lines) is very difficult. With a computer, such calculations 
are straightforward. Suppose you know a point, (x, y, z), on a flux line. 
(Any point (x, y, z) lies on some flux line. So choosing an arbitrary starting 
point defines a particular flux line which will be followed.) Suppose further 
that you want to take a step of size As along the flux line. The first 
defining property for flux lines (B is parallel to the line) implies that 


> 
Ax = As B,/|B| 
—> 
Ay = As BY /|BI 
— 
Az = As B,/|B| 


These equations say that As = x + AY + Az is parallel to B = By + By + 
B,. At the point (x, y, z) you calculate B: and then calculate Ax, Ay and 
Az. The next point along the flux line is (x + Ax, y + Ay, z + Az). Since 
you now know another point on the flux tine, you can repeat the process. 
In this way the computer is used to walk along the flux line in small steps, 


As. 





Flux lines never end, so the process continues until we return to the 
starting point. A block diagram of this strategy looks like this: 


INPUT 


Currents and Their Positions 


: > 
Choose a Point, r, 
on the Fiux Line 


Calculate the Total 
> 


Field, B, at r 


Calculate Ax = As B,/1BI 
> 
Az = As B,/|B| 


Caiculate new r 
= (x + Ax, y + Ay, z + Az) 





The fact that flux lines never end is a consequence of the fact that the 
next flux through any region must be zero (¢ B - dA = 0). As much flux 
must enter as leaves; hence, no flux line can end. (The corresponding 
statement for electrostatic fields is $ E + dA =4mkoq. This implies that 
electrostatic field lines end or start only on charges, which may, if 
necessary, be at infinity.) 


Notice that this iterative procedure for following flux lines will break 
—_> 
down if |B] = 0 at any point on the line. 
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FLUX LINES FOR A SET OF LONG STRAIGHT WIRES 


Let’s consider a simple case. 


Example 1. Map the magnetic induction pattern due to four long straight 
wires parallel to the z-axis. The currents and the points (x, y) at which the 
wires cut the x-y plane are: +1 at (+1, +1); -1 at (-1, +1); -1 at (-1, +1); 
and +1 at (+1, -10). 


Schematically the situation in the x-y plane is: 


Y 


QO | © 


(0,0) 


O | @ 


This is a relatively simple current distribution, but analytical methods 
cannot calculate the flux lines. A plus one current means a current of one 
ampere (in mks units) which is coming out of the page (by the right-hand 
rule). 


By symmetry, the magnetic field, B. has no z-component. All planes 
parallel to the x-y plane are equivalent, so we need only consider the x-y 
plane itself. A program implementing the strategy discussed above might 
look like the following: 


166 DIM XC4)3,YC0 413,104) 


118 LET N@s4 Currents and 


12@ FOR J=! TO N@ Geometry for 
138 READ X€J1,YC€CJUI],1CJd) long, straight 
140 NEXT J wires 


15¢ DATA loa les lewmlelewlewlerilsnmlels-ls } 

160 LET K@=16-E=07- Units 

172 LET D=.@5 — Stepsize, As 

18@ PRINT "(XsY) STARTING PT. ON FIELD LINE?" 
198 INPUT X5.Y5 

268 LET X@=x5-x 

21@ LET Y@=Y5-y 

220 LET X1=0-Ax 

230 LET Yl=@-Ay 

248 LET S@=.5 -— Print every S@ along line 

258 LET S=2 — Distance, S, along line 

262 PRINT xX", ¢y" 

278 LET Bl=a 

282 LET B2sg 

298 FOR J=1 TO NO 

300 LET RO=(XO+X1/2-X0d])*(XO+X1/2"KX0J1) i ie 
310 LET R2=R2+(YO+¥ 1/2-YCu}*CY OF 1/2-YCN)))/"! | Calculate B 
328 LET B=2eKQ«xICJI/R2_ 81/7 ihalt-stepped) 
336 LET BI=BI-CY@+Y1/2-Y(J])*B-8, 

340 LET B2=B2+(XO+X1/2-XC5) )*B- By, 

35@ NEXT J : 

368 LET B@=S@R(B1*B1+B2*Be) - (8! 

37G LET X1=D*B1/B@- Ax 

380 LET XO=X@+x 1—Newx 

390 LET Y1=D*B2/BB -— Ay 

400 LET Y@=2YO+¥ 1- New y 

416 LET $=S¢D 

426 IF S<S@ THEN 458 | print out 


436 LET S=@ Lines 
44@ PRINT XO,YO 
450 IF ABS(X@)+ABS(Y@)>1@8 THEN 47@ - Far off page oe 
460 IF ABS(CX@-X5)+ABSC(YB-Y5)>.9*"D THEN 279 - Back to | End of Line 
47@ PRINT X@,YQ@ — Last point on fine Start 

hs “#@ PRINT 
“9@ GOTO 188 - Return for new line 


48 END 
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The calculation of the magnetic field, Bir), in this program uses a point 
(approximately) one-half a step ahead of the present point on the line. A 
discussion of convergence and the half-step method are found in Appendix 
B. Arun of this program looks like this: 


RUN 

BMAP 

(XsY) STARTING PT- ON FIELD LINE? 
208408 

x Y 
162319 648183 
1- 49682 1-@6083 
128171 1- 38555 
©7439 17 945335 
792453 °804443 
(XeY) STARTING PT. ON FIELD LINE? 
725,90 

x Y 

° 382484 -0998777 
+ 49 69 48 ~1-4741 
798886 ~12¢87544 
1- 22665 -2- 14756 
1. 68888 -2.27552 
2e 18638 2026245 
2+ 66452 ~2e 12679 
3e 89 446 ~1+8678 
3¢ 45537 -1.52326 
3-7 323 -1419812 
30997 67 ~e 151479 
30977 37 ¢ 347 226 
3085484 °§31859 
3 63494 1.27863 
33226 1e 66852 
2093244 1297947 
2. 48 867 2- 19889 
1-99 878 20 28 379 
1.49415 20 24244 
1683185 2 657 38 
« 656388 1673115 
« 428 638 t~ 28974 
© 388341 794198 
« 471962 301455 
» 5860743 3e28055E-83 


(XeY) STARTING PT- ON FIELD LINE? 
? 


Practice Exercise: Modify this program to follow flux lines due to two 
wires both with | = +7 with (x,y) position (-.5,0) and (+.5,0). 


Terminal plotting is sometimes useful to display flux lines. If you have a 
plotter system or a plotting terminal, you may want to use that. Terminal 
plotting can be performed as discussed in Appendix A. The necessary 
changes in the flux line program BMAP on page 24 (and then a run) to do 
terminal plotting might look like: 


GET- BMAP 


10 READ X8,X9,Y8,Y9,Z28 
28 DATA -5,5,75, 551 

38 GOSUB 982398 

285 GOSUB 91862 

258 

430 GOSUB 9188 

45@ GOSUB 9162 

ase GOSUB 9200 
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RUN 
BMAP 


(XeY) STARTING PTe ON FIELD LINE? 

2eTan oT 

MAX Y= 5 
XXXXXXKXXXKAKKKXKKKKKKKXAKXKRARKAKAKKXKAKXAXXKXKKKXKXXXAXXXAXK 


Sa KR HK RK KKK KKH 


XXXXKXAKKAKRKKAXXKKXKKKXKKAXKAXKAXKKXAAKAKKXKXXX AX KXXKKK KKK 
MIN Y#-5 MIN X=<5 MAX X= 5 


(XsY>) STARTING PTe ON FIELD LINE? 

72558 

MAX Y= 5S 
XXXXXKXXXXXXXXXXXKXXXXKKXKXKXKXXXXKXKXKXXKXKKXX KK KXKKK KKK XKK 


SOR R RK KKK KKK RK 


XXKXAKKXKKXAXXKXKXK KKK XX KK KKX KKK AXXKK XX AXKAXKKXKXXKKXKK KKK KKK 
MIN Y3-S MIN X2=5 MAX X= 5 


(XsY) STARTING PTe ON FIELD LINE? 
? 


Practice Exercise: By hand, plot the positions of the line currents on the 
terminal plots. 


The flux lines bend around and close on themselves. The field is higher 
halfway between the oppositely-directed currents and lower halfway be- 
tween the like-directed currents. The field is quite uniform (that is, the 
flux lines are nearly parallel over a fair-sized region) at the center of the 
pattern. 


THE DENSITY OF FLUX LINES 


If the flux lines are started correctly somewhere tn the pattern, the density 
of the lines will be correct everywhere. That is, if the second defining 
property is obeyed somewhere by all the flux lines, it will be obeyed 
everywhere. The difficulty is that there is no easy way to make the density 
correct anywhere. (This is unlike the case for electrostatic field lines in 
Chapter One. For electrostatic field lines you can start the lines correctly 
near the charges.) 


You might be interested in thinking about this difficulty. Can you see a 
way to start the lines with the correct density? How about calculating the 
flux through small regions across, say, the line from (- 1,0) to (+1,0) in our 
example? You might then start the correct number of lines in each region. 


FIELD LINES FOR COMPLICATED CURRENT 
DISTRIBUTIONS 


The same algorithmic procedure to follow flux lines can be used for more 
complicated current distributions. You can use the Biot-Savart law quoted 
at the beginning of this chapter. By breaking up any current distribution 
into little pieces, lag, using the Biot-Savart law to calculate the dBs for 
each piece, and then adding the dBs to find the total magnetic field, B. 
you can calculate the magnetic field anywhere in space for any current 
distribution. Using the resulting values of B = (By, By, B,), you can follow 
the flux lines for any current distribution by making As parallel to B as 
above. 


Let's find the field pattern for a Helmholtz pair of circular coils. A 
Helmholtz pair has two parallel, circular coils carrying equal currents and 
having a separation equal to the coils’ (common) radius. 


Example 2. Map the magnetic field pattern in the x-y plane due to a pair 
of Helmholtz coils of radius one carrying a current of +1A. Both coils are 
parallel to the y-z plane. 
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> 
A program which calculates the magnetic field, B, using the Biot-Savart 
law and follows flux lines algorithmically might look like the following: 


COILMP 

1@@ DIM ClC26),SC€20) 

11@ LET N@=e2 

128 FOR J=1 TO N@ 

138 READ XCJI,RC0J},103) 

14@ NEXT J 

152 DATA ° 5, le lone Se ls? 

160 LET P2=6. 28318 

172 LET N9=16 

186 FOR Il=1] TONS 

19@ LET T=Pe*x C1l-.5)/N9 | Store 

206 LET CLII=COSCT) sines and 
21@ LET SCII=SINCT) cosines 
228 NEXT I 

230 LET K@=1].E~-@7 

246 LET D=.1 

25@ PRINT "CXsY) STARTING PT. ON FIELD LINE?” 
26@ INPUT X5,Y5 

278 LET X@=X5 

286 LET Y@=Y5 

292 LET X1=¢g 

306 LET Y1=2 

310 LET S@=.5 

320 LET S=6 

330 PRINT '%'', ''Y"* 

34@ LET BIl=@ 

350 LET. Be=@ 

360 LET B32@ 

37@ FOR J=1 TO N@ -—- step through coils 
38@ FOR 121 TO N9 — step around each coil 
3928 LET Li=@ 

A@@ LET Les-PeeRCJI*SCLI/N9 
41@ LET L3=PeeRCJI*CCLII/N9 

426 LET L7=xXtJ} 

43@ LET L8=RCJ1*CO1} 

440 LET L9O=PCJI*SC1) 

45@ LET X62X%@+X1/2-L7 

460 LET Y6=Y@4Y1/2-L8 

47@ LET Z6=-L9 

48@ LET R6=SQRCX6*X 6+Y 64*Y 6+Z 6*Z 6) 
490 LET R6=R6*«R6*RE 

502 LET Cl=L2*Z6-L3*Y6 

518 LET C2=L3*xX 6<-L1*Z os oe 
520 LET C3=L1*Y6-Le*xx6 

530 LET BI=Bl+K@*1CJ53*C1/R6- 8, 
540 LET Bez=Be+KO@xI CJ} *C2/R6-8B 
208: LET B3*B3+K Gel CJI*C3/R6- 8B, 
56@ NEXT I 

$7? NEXT J 

58@ LET B@=SQR(B1*B1+B2*B2+B3*E3) 
59@ LET X1!=D*«5)/B@ 

600 LET X@=X2+X1 

610 LET Y1=D*Be/Re@ 

620 LET Y@=YOQ+Y 1 

632 LET S$=S+D 

64@ IF S<S@ THEN 676 

652 LET 522 

662 PRINT X@sY@ 

67@ IF ALS(CXO)+ABS(YG)>1@8 THEN 692 
6683 IF ABS(CX@-X5)+ABS(YO-Y5)>.9*D THEN 342 
692 PFINT X@sY@ 

70@ PRINT 

718 GOTO 25¢e@ 

726 END 





Practice Exercise: Annotate the program by written comments by each 
line. 





A run looks like this: 


(XsY) STARTING PT. 


(XsY) STARTING PT. 


? 


RUN 
COILMP 
(X,Y) STARTING PT. ON FIELD LINE? 
78s 1 
x 
«469244 857836 
»789 897 1.18179 
« 615597 1.63534 
* 189641 1286649 
m2 312249 1-832@18 
~»7820194 1.52856 
--757882 16845849 
0 341947 © 869398 
~6-5458 6E-22 985284 


ON FIELD LINE? 


Qe 4s 0D 
x 
© 7 687 49 1613971 
¢ 596547 1.59283 
152524 1.885904 
--« 33987 11-7566 
-~o711667 1-43789 
~- 684021 °965799 
~+ 2269359 °©939789 
235256 932256 
« 325676 «889844 
€XsY) STARTING PT. ON FIELD LINE? 
27025 
x 
« 4997 65 «4973 
°97 69 « 635858 
be 37 464 °935883 
le 686935 1232485 
1692196 1.76524 
2- 08598 2e 237 29 
2- 17837 2272884 
221968 3¢ 22724 
2+ 13798 3272326 
1-998 62 4.289284 
1-77 646 4-65 
1.47184 584543 
1.89804 S- 3667 
« 644855 55-5921 
© 156582 5.69515 
=e 341449 5- 6788 
-e«816714 525197 
~it-e 24889 5225712 
-1-86911 4«488e22 
-2-@6825 4-e@269 
~2-16923 32539 49 
~2-e 19883 3° B84ass 
~2+ 15247 254345 
~22-83318 2e85834 
“1684265 1-59 658 
-1-57985 1217213 
-1+237 086 «809351 
~--826264 » 56855 
-.313676 «49 2488 
-1-38@63E-@2 «580097 


ON FIELD LINE? 
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And, finally, a terminal plot of the program, COILMP on page 28, is 


shown below: 


GET-COILMP 


1@ READ X8sX9,Y8,Y9s28 


20 DATA -5,.5,-5,5,28 


38 GOSUB 9882 


265 IF X5=999 THEN 715 
266 LET 28=26+1 


2685 GOSUB 918@ 
338 

668 GOSUB 93100 
69@ JOSUB 9188 
715 GOSUB 9288 


APP-TTYPLO 


RUN 
COILMP 


(XeY) STARTING 
282 1 


(XeY) STARTING 
Pete eD 


(XeY) STARTING 
72S 


(XsY) STARTING 
799938 
MAX Y= § 


PT.» ON 


PTe ON 


PTe ON 


PTe ON 


FIELD LINE? 


FIELD LINE? 


FIELD LINE? 


FIELD LINE? 


XXXXXXXXXKAXKRKRAAKXKAKXXKXAXKKXAKXKKXXKXKXXKXXKKK XXX KX AXKKKXKKXX 


Se K RK KKK KR RK KKK 


3 


3 


3 


3 
3 
3 
3 
3 
3 
1 
2° Be 3 
2 1 3 
l2l2ie2 2 3 
3 3 3 
3 3 


KAXKXXXXKKXXKRKXKXXAKAK XXX XKXKAKXXX KXAXKAKKKKAKXAKXKXXXKKXKXKXKAXA 
MIN Y#-5 MIN X#=5 


Again, the strategy used in this program is entirely general. The magnetic 


MAX X= 5 


>> 
field, B(r), for any current distribution can be calculated from the Biot- 


Savart Law: sets of short wire segments; open-wound, short solenoids; sets 
of square coils; anything at all. A number of possible cases are treated in 


the exercises at the end of the chapter. 


CONCLUSION 


This chapter has introduced a general method to map magnetic field 
patterns. Both simple sets of long straight wires were discussed and a 
completely general procedure based on the Biot-Savart law was presented. 
By running several situations you can get a good intuitive feel for the 


behavior of magnetic fields. 





CHAPTER TWO EXERCISES 


Consider two long, straight, current-carrying wires parallel to the 
z-axis. The currents and the points (x,y) at which the wires cut the 
x-y plane are: +1 at (1,0); +1 at (-1,0). 


a) Determine the magnetic field along the x-axis from -3 to 3. 


b) In the region .5 <x <1.5,-.5<y <5, calculate six magnetic 
flux lines such that the number of flux lines through a given 
area on the x-axis is proportional to the flux in that region. 
(Since the field is independent of z, the number of flux lines 
through a Ax should be proportional to the average magnetic 
field on that region.) 


Cc) In the region -3 <x <3,-3<y <3, calculate six magnetic 
flux lines such that the number of flux lines through a given 
region is proportional to the field in that region. 


d) Inthe region close to each wire, how does the field compare to 
that of a single current-carrying wire? In the region far away 
from all the wires, how does the field compare to that of a 
single current-carrying wire? 


Consider two long current-carrying wires parallel to the z-axis. The 
currents and the points (x,y) at which the wires cut the x-y plane 
are: +1 at (1,0); - 1 at (- 1,0). 


a) Determine the magnetic field along the x-axis from - 3 to 3. 


b) In the region .5 <x <1.5,-.5 <y <5, calculate six magnetic 
flux lines such that the number of flux lines through a given 
area on the x-axis is proportional to the flux in that region. 
(Since the field is constant in the z direction the number of 
flux lines through a Ax should be proportional to the average 
magnetic field on that region.) 


c) In the region -3 <x <3,-3<y <3, calculate six magnetic 
flux lines such that the number of flux lines through a given 
region is proportional to the field in that region. 


d) Inthe region close to each wire, how does the field compare to 
that of a single current-carrying wire? In the region far away 
from all the wires, how does the field compare to a single 
current-carrying wire? 


Consider the simple case of a single loop of current-carrying wire. 
Say the loop lies in the y-z plane centered at the origin with a radius 
of one and carries a current of one. 


a) Using the Biot-Savart law and the computer integration 
method, calculate the magnetic field at several points along the 
x-axis. (x = 0, 1, 5, 10). 
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b) It is possible to calculate these axial fields analytically. Do so 
and compare the analytical results with the ones obtained from 
the computer method. 


c) Evaluate the magnetic field along the y-axis out to the edge of 
the coil, 


d) Calculate the flux line pattern in the y-z plane through the 
coll, 


Consider the case of four current-carrying wires parallel to the z-axis. 
The currents and the points (x,y) at which the wires cut the x-y 
plane are: +1 at (1.5,1); -1 at (-1.5,1}); +1 at (.5,- 1); and -1 at 
(- .5,- 1). 


a) Sketch the flux lines close to each wire. Now sketch the whole 
flux line distribution. 


b) Calculate the flux line pattern for this distribution. 


Consider the case of four current-carrying wires parallel to the z-axis. 
The currents and the points (x,y) at which the wires cut the x-y 
plane are +3 at (2,0); -3 at (-2,0); -1 at (1,0) and +1 at (-1,0). 


a) Sketch the flux lines close to each wire. Now sketch the whole 
flux line distribution. 


b) Calculate the flux line pattern for this distribution. 


Consider a square Helmholtz pair of coils. Let the corners of the 
square loops be defined by the (x,y,z) points: (+.5,1,1), (.5,1,-1), 
(.5,-1,-1),  (.5,-1,1) and (-.5,1,1),  (-.5,1,-1),  (-.5,-1,-1), 
(-.5,-1,1). 


a) Calculate the flux line pattern for this distribution. Since it is 
not rotationally symmetric about the x-axis, calculate it in the 
x-Z plane and then in the plane defined by the point (0,1,1) 
and the z-axis {i.e., the plane at 45° to both the x-z and y-z 
planes). 


b) Compare the axial field to that of the true (circular) 
Helmholtz pair. Compare both the magnitude of the field and 
the uniformity of the field in the center of each pair. 


Consider an infinite current-carrying cylinder of radius one. Model 
the cylinder as sixteen infinitely long wires parallel to the z-axis. 


a) Calculate the field distribution for this pattern. 


b) How well does this compare to the Biot-Savart law inside the 
cylinder? Outside the cylinder how does it compare in magni- 
tude and shape to a single wire at the origin with a current of 
sixteen? 





10. 


The field for an electric dipole was given in Chapter 1, Exercise 10. 
How well does the magnetic field of a simple loop in the x-y plane of 
radius 1 compare to the dipole equation when you are far away from 
the loop? Use the distance of ten for your comparison. What dipole 
moment, p, Is necessary to make the magnitudes agree? How does 
this dipole moment vary with the current in the loop? 


Flux lines are not the only method of studying magnetic fields. 
Write a program similar to the equipotential character mapping 
program to print like characters for regions with the magnitude of 
[B| equal. This method could be useful for determining regions of 
constant field. 


a) Use your program to study the field pattern in Exercise 7. 

b) Use your program to study the field pattern in Exercise 2. 

Consider a solenoid of radius 1 along the x-axis five units long 

centered at the origin. (Approximate the solenoid by five current- 

carrying loops). 

a) Calculate the field line distribution. 

b) — Using the field line mapping program or the field magnitude 
program of Exercise 9, determine what current should be run 


through the end coils to make the field inside the solenoid 
more uniform, 
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CHAPTER THREE: THE MOTION OF 
CHARGED PARTICLES IN ELECTRO - 
STATIC AND MAGNETOSTATIC FIELDS 


INTRODUCTION 


Many students believe that charged particles move along field lines. Except 
in very rare instances, charged particles do not move along field lines. 
Although the direction of the field line at any point defines the direction 
of the force at that point, the force only changes the velocity not the 
position itself. Once it is moving, the charged particle will, in general, cross 
field lines because of inertia. (If the particle moves through a very viscous 
medium, you can force the particle to stay on a field line.) 


This chapter introduces a way to calculate the trajectories of charged 
particles in any (combination of) electrostatic and magnetostatic fields. 
Those of you who have seen p= ma algorithmic solutions with a computer 
before (probably in mechanics) will recognize the discussion in this 
chapter. The chapter is based on an , algorithmic solution for particle 
trajectories when the forces are due to E and B. 


The forces on charged particles in electrostatic and magnetostatic fields are 
discussed in physics textbooks. _the force on a charged particle in an 
electrostatic field, E, is just F(r) = qE(r), where q is the value of the 
charge. The force on a charged particle in a magnetic field (or, more cor- 
rectly, magnetic induction), B. is the Lorentz force, F(z) = q(v(r) x Bir). In 
combined electrostatic and magnetostatic fields, the force is F = q(E+vxB). 


AN ALGORITHMIC SOLUTION TO NEWTON’S SECOND 
LAW 


E = ma is Newton's Second Law and is the basis of classical calculations of 
particle trajectories. In most physical situations, you know the force, Fir). 
Hence, you know the acceleration a(r) = F(r)/m. Given a way to find the 
acceleration at any point in space, the computer can easily be programmed 
to find the trajectory of any particle. 


Suppose you know a point, f= (x Y,z), on the trajectory of a particle and 
also Bnew the particle’s velocity, v= (Vy,V yr Vv ;) at that point. (You usually 
know rf and ¥ at the initial point of the trajectory.) Suppose you want to 
find the position a short time, At, later. From the acceleration, you can 
find the change in the velocity of the particle in the time, At: Av = at. 
Then, the new velocity after the time, At, is 


‘ = 
Vie =Voig + Av (8) 


From the velocity, you can find the change in position r during the time 
At, since 


r=vAt (9) 


36 O ELECTRIC AND MAGNETIC FileLOS 


Hence, the next position on the particle’s trajectory is 


t soe = Fold + vAt (10) 


The new time is they = tojg + At. 


Now that you know another point on the trajectory, the process can be 
repeated over and over. The result is the trajectory of the particle as it 
moves through the force, F(r). We will discuss two cases: the first is the 
motion of charged particles in electrostatic fields, and the second is the 
motion of charged particles in combinations of electrostatic and magneto- 
static fields. The second case is somewhat more complicated because the 
force on a charged particle in a magnetic field is velocity-dependent. 


THE MOTION OF A CHARGED PARTICLE IN AN ELEC- 
TROSTATIC FIELD, E 


The force on a charged particle in an electrostatic field, c. is just F = gE. 
You can calculate the electrostatic field, E for any distribution of point 
source charges and line source charges. The field for a point source charge 
iS E = Ko O/r? (t/r) where Q is the value of the source charge and r is the 


distance from the source charge, r = irl. The field for a line source charge ts 
-> 
E = kgA/r (r/r) where A is the charge per unit length of the line charge and 


r is the perpendicular distance from the line charge. ky defines the units; 
ko = 9 x 10? in mks units and kg = 1 for cgs units. Any general source 
charge distribution can be treated as a set of point or line charges. 


Example 1. Find the trajectory of a positron (a particle with the mass of 
an electron but opposite (positive) charge) moving under the influence of a 
uniform E field parallel to the y-axis. Start the trajectory at the origin with 
velocity (4,0,0). 


One program which implements the algorithmic solution discussed above is 
shown below: 


MT ONE 


1@@ READ E1lsE2,E3 os 

118 DATA @51,8 | see & rea 

120 PRINT “INITIAL XsY¥sZs VX, W.VZ & FINAL eS Initial 7, Vand 
138 INPUT XsYsZs Vis V2. V34T2 final time 
142 LET De.@5-At 

15@ LET Q@=1-—aq/m 

168 LET Tx@ 

178 LET T9238 

18@ REM COMPUTE FIELDS, ACCELS., VELes, POSITION, TIME 
198 LET Al®Q@*E1-— 4, 

20@ LET A2=Q0*E2-a4 

21@ LET A3=Q0*E3-a, 

226 IF T>@ THEN 276 

238 LET Vit=V1I4Al*"D/2 

248 LET V2=V2+A2«D/2 | initial half step 

258 LET V32V34A3*«D/2 

260 GOTO 38a 

e708 LET VIFVI+Al*D-vV, 

260 LET V2=Ve+AexD- Vv 

2986 LET V3=V3+A3*«D- Vv, 

306 LET X2X+V1*D_—NewXx 

318 LET Y2¥+V2*eD- New Y 

320 LET 222+V3*D—NewZ 

338 LET T=T+D—- Newt 

342 LET T9=T9+D 

358 IF ABSCT9-T9/18)/T9>.2@802081 THEN 382 | Print results 
362 LET TO2#Q every TO/10 
37@ PRINT T3Xs3Y3Z 

380 IF T<TQ THEN 18@ — Return for next point on trajectory 

398 PRINT 

4@Q GOTO 128 — Return for new trajectory 

418 END 
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The program uses a half-step method (the velocity is one-half a time step, 
At, away from the position). This half-step method increases the con- 
vergence of the procedure and is discussed further in Appendix B. The 
units used in the program are normalized (because q/m is set to 1) so that 
fields are measured in units of the charge-to-mass ratio of the particle. For 
example, since the particle is a positron, q/m = +1.759 x 10!! in mks 
units, so a field with |E| = 1 actually has |E| = 5.685 x 10'? V/m. 


A run of this program looks like the following: 


RUN 
MTONE 


INITIAL XsYsZsVXsVY,VZ & FINAL T? 
2056s 8s 458,28 18 

le 4-6 2-5 B 

2e 7699999 Be G 

Ge 126 45 8 

4.@09001 166 3. @ 

5eOGG81 206 12-5 @B 
6.08801 24.0001 18. @ 
7+O82@882 28-80281 24-5 @ 
8-@89882 32-9801 32-.A881 B 
9-O@08Al 36- 48.580! @B 
9-99999 49. 52-0081 2 


INITIAL XsYsZs VX,VY2VZ & FINAL T? 
? 


You can plot the trajectory of the particle using MTONE and the terminal 
plotting discussed in Appendix A. The required changes in the program 
and a run of the modified program are shown below: 


GET-MTONE 


18 READ X8sX9,Y8,Y9,ZG0 
28 DATA @.58+ 8,585 1 

38 GOSUB 9888 

132 LET X@=X 

134 LET Y@zyY 

136 GOSUB 9162 

378 LET X@2X 

372 LET Yaryv 

374 GOSUB 9188 

398 GOSUB 9288 


APP-TTYPLO 


RUN 
MTONE 


INITIAL XsYsZoVKseV¥.VZ & FINAL T? 

78s Os Gs 4o Bs Bs 1B 

MAX Y= 5@ 
XXXXXXXXKXXXXKXKXKKXKXKKKKKKXKKXKKK XXKKKXKKXK XK KK XK KK XX KK 


1 


See CHR RK KKK RRR KKK KKK 


hat 


I 
AXXAXXXXXXXXXAXXXXKXXXXKAXKXXAXXXKXX KXAXKXKKXKXXKXXKAXKKKXXX KXXKXKK 
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MIN Y= 6 MIN X= @ MAX X# 58 
INITIAL XsYsZsVXsVY,VZ & FINAL T? 
f 


For a uniform electrostatic field, the trajectory is a parabola (as can be 
derived analytically). Consequently, the motion is just like that in a 
uniform gravitational field. Thus, this particular example can be used to 
check the numerical method. This check is pursued in a problem at the 
end of the chapter. 


Practice Exercise: Modify the program, MTONE, so that it calculates the 
electrostatic field inside the integration loop. Use the electrostatic field 
due to a point charge at the origin. 


A terminal plot of several orbits in the field of a -1 charge at the origin 
looks like: 


16 READ X8sX9sY8sY92Z0 
28 DATA 20 39°225, 2-5 
36 GOSUB 9898 

182 

118 

131 IF X2999 THEN 4¢@5 
132 LET Z202Z8+1 

133 LET X@2X 

134 LET Y@sv 

13S LET R3eCCX1+X1/2)t Se¢CVY 14V1/2)12)¢91-5 
196 LET Al?@ 

206 LET Ae=6 

218 LET A3=e@ 

37@ LET X@=sxX 

372 LET Y@=Y 

374 GOSUB 91286 

398 

485 GOSUB 9222 


APP-TTYPLO 


RUN 
MTONE 


INITIAL XsYsZs VX, VY,VZ & FINAL T? 

7-1,0. Gs Bs ts a7 

INITIAL XsYoZsVXeVYsVZ & FINAL T? 

2<1,65 8,05 1-25 Bs 15 

INITIAL XsYsZsVX,VYsVZ & FINAL T? 

79995 Bs Os Gs Os Os @ 

MAX Y# 2-5 
AKRXKKAXXAKKAKXXKXKXKXXKXXXXKKXXKXXKXXXXKXKXXXKXKXKXXKXKXXKXKXXK 


Y 
Y 

Y 

Y 

¥ 

Y 2 

Y 

Y 2 2 

Y 

Y 1 

Y i 

Y t . 2 
Y 

Y ! 

Y 

Y 2 2 
Y 

Y i 

Y 

Y ! 2 
Y 1 

Y ! 

Y 

Y 2 2 

Y 

Y 2 

Y 

Y 

Y 

Y 

Y 


XXXXKXKXAXXAKXAKARKXKXKARKKX KXKXKXXAXKXKKKXKAAKXX KK KK KKXKAXXKKXKKXKXKXKKX 
MIN Y3-2e5 MIN X==-2 MAX X= 3 


END 
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Practice Exercise: Plot the position of the source charge on the terminal 
plot. 


These trajectories can also be calculated analytically, although with some 
difficulty. The trajectories are exactly similar to the orbits of a satellite 
around the earth. 


Practice Exercise: Modify the program to find the trajectories of a posi- 
tron in the electrostatic field due to a +7 point charge at (+.5,0,0) anda-1 
point charge at (-.5,0,0). 


These trajectories cannot be derived analytically. Other source charge 
distributions are considered in the problems at the end of the chapter. 


THE MOTION OF CHARGED PARTICLES IN COMBINED 
ELECTROSTATIC AND MAGNETOSTATIC FIELDS 


The algorithmic solution for charged particle trajectories becomes more 
difficult when the force is velocity-dependent. The program discussed 
above for electrostatic fields calculates the velocity a half-At time step 
away from the position and the acceleration. Since the acceleration needs 
the value of the velocity at the acceleration’s own time, an estimate of this 
un-half-stepped velocity must be made. This estimated velocity can then 
be used to find the acceleration. 


Example 2. Find the trajectory of a positron starting at (-1,0,0) with 
velocity (0,1,0) in a uniform magnetic field of magnitude (q/m) parallel to 
the z-axis. 


The following listing is for a program which implements the algorithmic 
trajectory calculation for charges moving in combined electrostatic and 
magnetostatic fields. 


MTONEB 


168 READ E1,E2,£E3.,.B1., Be, B3 

118 DATA Gs Bs Gs Bs B» 1 

120 PRINT “ENITIAL XsYsZsVXsV¥s.VZ & FINAL T?" 
132 INPUT XsYs.Zs. Vile Ves V3 TS 

140 LET D=.85 

156 LET Q@=1 

168 LET Ts 

170 LET T92e@ 

1$0 LET Al=8 

19@ LET Ae=2 

282 LET A3=8 

21@ REM COMPUTE FIELDS, ACCELS+e, VELes POSITION, TIME 
22@ LET UI=V1I+Al*D/2 

230 LET UezV2+A2e*«D/e2 

240 LET U3=V34A3"D/2 

25@ LET Al=QO*(E1+¢( U2*B3-U3*Be) > 
260 LET A@™=Q0*( E2+(U3*4B1-U1*B3) > 
270 LET A3=Q0*(E34+(U1*Be-Ue2*Bl)) 
28@ IF T?@ THEN 336 

298 LET VI=V1+Al*D/2 

398 LET Ve=Ve+Ae*D/2 

31@ LET V3=V34+A3*"D/2 

328 GOTO 368 

338 LET Vi2eVi+4aAl*D 

348 LET Ve=V2+Ae*"D 

358 LET V3=V3+A3*D 

366 LET X=X+V1*D 

376 LET Y2¥+Vex«xD 

382 LET 23Z+V3*D 

3980 LET T=T+D 

4@@ LET T9=T94+D 

418 IF ABSCT9-T6/18)/T9>-A@GGG6Gi THEN 448 
426 LET T9=6 
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43@ PRINT TsX,YsZ 
448. IF T<T@ THEN 216 
458 PRINT 

468 GOTO 129 

47@ END 


Practice Exercise: Annotate the program by written comments next to 
each line. 


A run of this program for the initial conditions given is: 


RUN 
MTONEB 


INITIAL XsY¥sZ.VXsVY,VZ & FINAL T? 
2-1, 00 8s Bs 158,18 


Le -2539729 *841534 8 
ce 417183 ©9879 16 a 
Se 989278 ° 1379 48 g 
4.02001 « 64953 <0 759143 @ 
Se 80801 -- 28915 -0956982 Q 
6+ 88081 “962105 --2273279 8 
7+@86Be -~e7T49 406 + 662267 6 
8- 8882 © 153869 987741 @ 
9+ 82891 914821 « 403429 @ 
9099999 832466 72552571 @ 


INITIAL XsYsZsVXeVY,VZ & FINAL T? 


2 


The motion is a circle as is expected from the analytic solution. A sensitive 
test of the accuracy of this computer-based method is an examination of 
the magnitude of the velocity as the particle moves around the circle. The 
matter of accuracy is examined further in an exercise at the end of the 
chapter. 


Motion in a magnetic field often involves all three directions, x, y and z. 
Hence, the motion is often hard to plot on a two-dimensional piece of 


paper. 


Practice Exercise: Modify the program to compute trajectories when both 
a untform electric field in the y direction and a uniform magnetic field in 
the z direction are present. 


Such a field configuration is called a velocity selector. When the initial 
velocity of the charged particle has magnitude IEI/\B| and is perpendicular 
to both fields, the trajectory is a straight line. Thompson used this field 
arrangement in his early measurement of (e/m) for electrons. When the 
velocity is not perpendicular to both fields, the motion can be quite 
complicated. 


Practice Exercise: Modify the program to find the trajectories of a posi- 
tron when both the electrostatic field due to a point charge at the origin 
and a uniform magnetic field in the z direction are present. 


Analytical methods won’t produce these trajectories even though the 
changes necessary in the algorithmic program are minimal. Further physi- 
cal field situations are considered in the problems at the end of the 
chapter. 
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CONCLUSION 


This chapter introduced a way to calculate the trajectories of charged 
particles in electrostatic and magnetostatic fields. The method works for 
arbitrary fields and for any charged particle. The method is general but 
easy to understand. By running several examples you can gain a good 
understanding of the motions of charged particles in e & m fields. 


CHAPTER THREE EXERCISES 


1. Consider the motion of a positron in the electric field created by an 
infinitely long line charge. The charge is along the z-axis and has a 
magnitude of +1 coulomb/m. 


a) 


b) 


c} 


d) 


e} 


Determine the trajectory of a positron initially travelling in the 
Z direction starting at (1,0,0) with a velocity of one. 


Determine the trajectory of a positron initially travelling with 
a velocity of (5,0,0) and an initial position of (-40,7,0). 


Does the angle of deflection of the positron change with the 
different y values for its starting point in part (b)? Try y=.5, 


y=2, and y=5. 


What happens when the charge on the line is changed to the 


_ charge on a (negative) electron? 


Will the electric field ever be able to hold the positron (i.e., 
will the positron ever orbit the line charge)? 


2. Consider the motion of an electron in the electric field created by 
two infinitely long line charges. The charges are parallel to the z-axis, 
each with a line charge density of -1. The line charges intersect the 
x-y plane at the points (0,2) and (0,-2). 


a) 


b) 


c) 


d) 


Determine the trajectory of an electron initially travelling in 
the z direction starting at the origin with a velocity of one. 


Determine the trajectory of an electron initially travelling with 


' a velocity of (5,0,0) and an initial position of (-20,0,0). 


Does the angle of deflection of the positron change with the 
different y values for its starting point in part (b)? Try y=-17, 
y=-.5, y=.5, y=7. 


What happens when the charge on the line is changed to the 
charge on one (negative) electron? 
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Ke Consider the motion of a positron in the electric field created by 
two point charges. The charges and positions are: ~1 at (3,0,0) and 
-1 at (-3,0,0). 


a} What happens to the positron with the initial conditions: 
position (0,8,0) and velocity (.45,0,0)? 


b) What happens to the positron with the initial conditions: 
position (4,0,0) and velocity (0,1,0)? 


c) What happens to the positron with the initial conditions: 
position (3,0,3) and velocity (0,.5,0)? 


d) Try other trajectories that you think might be interesting. 


4, Consider the motion of a positron in the electric field created by a 
dipole charge distribution. The charges and positions are: -1 at 
(3,0,0) and +1 at (-3,0,0). 


a) What happens to the positron with the initial conditions: 
position (0,8,0) and velocity (.45,0,0)? 


b) What happens to the positron with the initial conditions: 
position (-4,0,0) and velocity (0,1,0)? 


c)} What happens to the positron with the initial conditions: 
position (4,0,0) and velocity (0,1,0)? 


d) What happens to the positron with the initial conditions: 
position (3,0,3) and velocity (0,.5,0)? 


e) Try other trajectories you think might be interesting. 


2. Constant magnetic fields are sometimes used as simple electron beam 
collimators. Suppose we have a constant magnetic field in the z 
direction of magnitude (e/m). 


a) What happens to an electron with initial position at the origin 
and initial velocity (0,0,17)? 


b) What happens to an electron with initial position at the origin 
and initial velocity (0,17,7). 


c) What happens to an electron with initial position at the origin 
and initial velocity (0,1,0). 


d) How does this field act as a simple collimator? 


6. Consider the diverging field By = 0, B,, = 1/x?, and B, = 1/x?. 
xX y Zz 


a) What happens to a positron with initial position of (20,0,0) 
and initial velocity of (-1,0,0). This is often called a magnetic 
mirror because the field pattern reflects charged particles. Two 
such mirrors contain charged particles between them and are 
often called magnetic bottles. 
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b) | What happens to the velocity in the x direction? 
c) /s energy conserved? 


Consider the combined electric and magnetic field problem discussed 
in the chapter (the velocity selector with a uniform E = q/m in the y 
direction and a uniform B = q/m in the z direction). 


a) Verify that those particles tha t are perpendicular to both ra 
> > 
and B and have a velocity of \|E|AB| pass through this field as if 
the field did not exist. 


-> -> 
b) In what direction do the particles move if \E|/AB| is greater 
than or less than the velocity? 


c) In what direction do the particles move if the x-component of 
> —> 
velocity is |E\/\B| but the electron has some velocity in the y 
direction? 


d) In what direction do the particles move if the x-component of 
—> — 
velocity is |E|A8| but the electron has some velocity in the z 
direction? 


-> — 
e) Determine analytically that V, = \E\/A8| for the velocity 
selector. 


Consider a uniform magnetic field of .1 in the z direction and a 
point charge of +1 at the origin. 


a} lf an electron has an initial position of (1,0,0) and velocity of 
(0,1,0), what does its trajectory look like? 


b) — If. an electron has an initial position of (1,0,0) and velocity of 
(0,0,1}, what does its trajectory look like? 


The uniform electric field problem is similar to the acceleration of a 
body under the influence of gravity. 


a) Write down the analytical expression in component form for 
the position of a particle as a function of time fincluding its 
initial velocity). 


b) Compare the results of the program MTONE to the analytical 
expression at the times t = 2, 5, and 10. 


The solution to the problem of a positron moving in a constant 
magnetic field is a good problem to check the iterative solution 
because there is an analytical solution. Start the charged particle 
with initial velocity perpendicular to the magnetic field. 


a} Determine (analytically) the radius of the (circular) orbit asa 
function of the charge, mass, magnetic field, and velocity. 
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b) 


c) 


d) 


Determine the period of the motion as a function of the 
magnetic field, charge, and mass. (The frequency associated 
with this period is called the cyclotron frequency.) 


Compare your results with that obtained by the algorithmic 
method. 


-> 
The analytical solution has |V| = constant. Check your iter- 
ative solution by examining how constant |V| remains as time 
passes. 





CHAPTER FOUR: GAUSS'S, 
AMPERE’S, AND FARADAY’S LAWS 
FOR STATIC FIELDS 


INTRODUCTION 


Gauss’s, Ampere’s and Faraday’s Laws represent the original forms of 
Maxwell's four equations. Aside from the term Maxwell added to Am- 
pere’s Law (which made the set of equations agree with conservation of 
charge and also predict electromagnetic waves), these four equations 
represent the basis of all classical electricity and magnetism. 


This chapter shows you a way to use the computer to calculate the 
integrals in Gauss’s Laws for E and for B, Ampere’s Law and Faraday’s 
Law. You can then deal with questions such as “Can this given field be an 
electrostatic field? If so, where are the charges and what values of charge 
are present?” or ‘‘Can this given field be a magnetostatic field, and, if so, 
where and of what size are the currents?” (Again we will use the term 
magnetostatic field for B even though magnetostatic induction is strictly 
correct.) In the forms we shall! use, the laws can be written as follows: 


—> > > 

Gauss’s Law for E fE*dA = 4rkoq (11) 
> > > 

Gauss’s Law for B $#B°dA=0 (12) 
> 

Ampere’s Law §B de = ky | (13) 

Faraday’s Law fE+ ag =Q (14) 


(In the general case Faraday’s Law says that §E-de = -do/dt, but we are 
dealing with static cases. Currents are constant in time so the flux, ® does 
not change.) The symbol ¢ means the closed integral: in Gauss’s Laws the 
closed integral is taken over some surface completely enclosing some 
volume of space; in Ampere’s and Faraday’s Laws the closed integral is 
taken along some contour completely enclosing a surface in space. ko 
defines the units, and gq and | are the net enclosed charge and current, 
respectively. 


These laws give a way to examine any field to see if it can represent an 
electrostatic field E (in which case $E*d2 must be zero everywhere and 
fE*dA can tell you the values of the charges present) or a magnetostatic 
field B (in which case $B*dA must be zero everywhere and §B-de can tell 
you the values of the currents present). These integral forms of the four 
basic laws also determine the four differential forms for Maxwell's equa- 
tions which use vector derivatives (the divergence and the curl). 


The chapter starts by presenting a number of pictures representing possible 
vector fields. These pictures will be examined to see which might represent 
electrostatic or magnetostatic fields. The chapter will then move on to 
show you how the computer can be used to deal more quantitatively with 
fields using the four basic laws in their integral forms. The most general 
question of how to deal with fields if you are given only a set of measured 
values throughout some region of space is considered in the exercises at 
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the end of the chapter. The discussion in the chapter will be given in terms 
of E. and B fields, but the methods are applicable to any vector field (such 
as the gravitational field or the velocity field of a flow pattern). 


PICTURES OF VECTOR FIELDS 


It is sometimes helpful to see a number of pictures of possible vector 
fields. Inspecting the pictures to determine whether or not each picture 
might represent an electrostatic or magnetostatic (or gravitational or 
whatever) field exercises your knowledge and comprehension of the basic 
laws. By inspection, you can often guess whether or not a particular 
picture could represent an E or B field; you can often guess where the 
source charges or source currents would lie. To be more quantitative, say 
to find the values of the charges or currents, you need to perform the 
necessary integrations. The computer helps in two ways: first, the com- 
puter can be used to produce the pictures of the fields, and second, the 
computer can perform the integrations. 


There are two common kinds of pictures which are used to represent 
vector fields for discussion. The first is a figure that shows the vector field 
represented as arrows at a number of points throughout a region. The 
second is a field line map such as those discussed in Chapters One and 
Two. We will concentrate on pictures of the first kind. Let us look at 
several examples: 





Figure (a) 
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Figure (d) 
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Figure (f) 





Figure (g) 
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Figure (i) 








Which of these figures could represent an electrostatic field? You need to 
identify those fields for which $E-de = 0 for all (closed) paths. Figure (a) 
is such a field; it appears possible that adding up E-Ag for small pieces 
around a contour made up of radial and circumferential pieces would 
result in zero. In principle, you would have to show that the integral was 
zero for any possible contour; in practice, if one type of contour has 
§E d2 = 0, it is likely that almost any contour will also have §E- de = Q. In 
the case of Figure (a), it seems clear that any contour which can be 
approximated by a Set of radial and circumferential pieces might obey 
Faraday’s Law (¢E*d& = 0}, and that is certainly a wide class of contours. 


On those figures you believe could be electrostatic fields, identify the 
positions and relative values of the source charges. In Figure (a) there 
seems to be a positive charge at the center of the figure (and an equal 
amount of negative charge at °° on which to terminate the field lines). 
When several charges are present, you can estimate the relative magnitudes 
and the signs of the charges by considering the field near each charge. 


Which of these figures could represent a magnetostatic field? This is a 
harder question because you really need to find fields for which §B- dA= 0 
and the figures are only two-dimensional. Let us agree that the physical 
situations represented in the figures have no dependence upon z, that is, 
that every plane parallel to the x-y plane is equivalent to the x-y plane. 
Then the x-y plane shown in the figures tells the whole story. Figure (b) 
> —> 

seems likely to have $B*dA = O for surfaces which are pieces of circular 
cylinders with zaxis axes connected by radial planes. Again you cannot 
prove that $B*dA= 0 for all possible surfaces without doing some integrals, 
but the result looks right. 


On those figures you believe could be magnetostatic fields, identify the 
positions and relative values of the source currents. In figure (b) there 
seems to be only one source current at the center of the picture. The 
current appears to come out of the page (by the right-hand rule), so the 
current is positive although its magnitude cannot be determined without 
more information. If there were several currents, their relative values could 
be estimated from the relative strengths of the fields nearby; the directions 
of the currents could be determined from the directions of the fields 
nearby. 


Now you can try the rest of the figures. It is not difficult for anyone to 
produce pictures such as those shown. You can do it yourself. 


> > 
Practice Exercise: Write a program to calculate \E| and the angle E makes 
with the x-axis for a number of points on the x-y plane. You may assume 
EYE y E, are given as equations. 
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SURFACE AND LINE INTEGRALS WITH THE COMPUTER 


Now let us see how to be more quantitative. It is easy to do the integrals 
involved in Gauss’s, Ampere’s and Faraday’s Laws on the computer. An 
exact integration would have to be analytical, but most field patterns 
cannot be integrated analytically. The method of integration we will 
present is approximate (although the method can be made accurate as 
desired) but it can integrate the necessary surface and line integrals for any 
vector field. 


To perform the line integral, you break the contour up into small seg- 
ments, Ag, evaluate the field, F, at the center of each segment, and form 
the dot-products Fe AX. Adding up the terms F°Ag for all the segments that 
form a closed loop yields an approximation to §F° dk for the v vector field 
F. Again, in the limit of the sum of infinitesimal segments, ds, the result 
would be exact. 


To perform the surface integral, you break the surface up into small areas, 
KA, evaluate the f field, E. at the center of each § small area, and form the 
dot-products FAA. Adding up the terms Fs KA for all the small areas 
forming a closed surface is then an approximation to the integral | $F: dA 
for the vector field e In the limit of the sum infinite simal areas, dA, is by 
definition the integral; for finite AA, the sum is an approximation. 


>— > 
Example 1. Evaluate $FedA for the vector electrostatic field, F, due to a 
line charge coinciding with the z-axis. Consider rectangular parallelepiped 
surfaces. 


A program to find this surface integral might look like the following: 


10@ PRINT "€XsY¥.2Z) FOR 2 | Set region 
118 INPUT X1,Y1]1,Z21,4X2sY2,Z22 

120 LET P@s@- Fiux 

136 LET NO@216 — #of segments in x, y andz 

140 LET K@=1- Units 

1820 LET Q=+ 1— Magnitude of charge 

168 LET XO=(X2-X1)/N@- Ax 

17@ LET Y@=¢(vy2-Y1)/N@—- dy 

188 LET Z26#(Z2-21)/NQ —- Az 

19@ LET A@=X@*Y@-— AA 

200 FOR X2X1+X6/2 TO X2-XO/2 STEP XB 
210 FOR Y=Y14Y0/2 TO Y2-Y@/2 STEP Y@ 
22@ LET F3RKO*«Q*Z2/( X*X+Y*Y+Z 24%Z 2) 
238 LET PO®PO+F3* AB 

240 LET F3OsKO*Q4Z 1/CX*X4+Y*Y4Z O*Z 2) 
258 LET PO=PO-F3*Ad 

260 NEXT Y 

27@ NEXT X 

260 LET AB=YO*ZG—AA 

290 FOR Y2Y1+Y@/2 TO Ye-YG/2 STEP Y@ 
380 FOR Z2*Z1+Z0/2 TO Z2-Z8/2 STEP ZG 
318 LET FlSKO*"Q*X2/(X2"X S+Y*XY4+Z*Z) 
326 LET P@=PO+F1* Ad Two Y-Z 
330 LET Fi=KO*xQeK1/ OX 1#X 1 4YRY4+Z*Z ) Planes 
3409 LET P@=Pe-F1* Ad 

356 NEXT Z 

366 NEXT Y 

378 LET AG*XO*ZB-AA 

38@ FOR X2X1+#xX@/2 TO X2@-XB/2 STEP XB 

39@ FOR Z=Z1+Z0/2 TO Z2-Z28/2 STEP ZO 

408 LET F2akKOx Q¥Y O/C X*X+Y S*Y 24+Z2%Z ) 

410 LET POs PO+F ee AS 

420 LET FOeKOeQRY 1 /CX#X+Y L#Y 14Z*Z) os 
43@0 LET P@2PG-Fe* ag 

44@ NEXT Z 

450 NEXT X 

466 PRINT "FLUX ='3 P@ 

4760 PRINT 

482 GOTO 188 

498 END 


Two X-Y 
Planes at 
Z?! and Z2 





The field, F = (F1,F2,F3), is the electrostatic field due to a line charge at 
the origin [koA/r(r/Irl)], SO $FdA should be zero unless the region 
enclosed by the rectangular parallelepiped encloses the z-axis. When the 
z-axis is enclosed, §F° dA should equal 47k gA(Z2-Z1). A run of this 
Gauss’s Law program looks like: 


RUN 
GAUSS 


(XsYsZ) FOR 2@ CORNERS? 
Flee dents ls lel 

FLUX = 1523683 

(XsY2Z)> FOR 2 CORNERS? 
Tha ly le2s2,2 

FLUX = 2.217896 

(Xs¥eZ) FOR 2 CORNERS? 
Pleo le ls 2s 203 

FLUX 2 2.355734 


(Xs¥sZ) FOR 2 CORNERS? 
? 


Practice Exercise: How much charge lies along the line charge between Z = 
O and z = 2? What is the linear charge density, }? Does this value agree 
with the equations in the program? 


Practice Exercise: Modify the program for line charges having \ = +1 and 
-1 which are parallel to the z-axis and cut the x-y plane at (+.5,0) and 
(-.5,0), respectively. Where should the surface integral be zero and where 
non-zero? 


Practice Exercise: Modify the program for a long, straight wire carrying 
* « e = = — * 

current | which coincides with the z-axis. (Remember that B is perpen- 

dicular to fr. }) Where should the surface integral be zero and where non- 

zero? 


Further examples are considered in the exercises at the end of the chapter. 
The field, F, can be given as equations or even as data. 


—_> 
Example 2. Evaluate §Fede for the vector magnetostatic field due to a 
long, straight wire carrying a current of 1 ampere and coinciding with the 
z-axis. Consider rectangular contours lying in the x-y plane. 
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A program to evaluate the line integral might look like: 


AMPFAR 


180 PRINT "*€X,Y) FOR 2 CORNERS?" 
118 INPUT X1L,Y1,X2sY2 

128 LET L@=zG 

130 LET N@=#5ie2 

148 LET K@#1-E-987 

158 LET 121 

160 LET X8=(X2-X1)N@ 

178 LET YO=(Y2-Y 1) /NB 

18@ FOR X*X1+X8/2 TO X2-X@/2 STEP X@ 
190 LET Fe-2eK Oel«Y 1 /CX*#X+Y 1*Y 1) 
220 LET LOeLO+raxa 

216 LET Fa-2aKQexl«YO/(X*"X+4+Y Quy 2) 
226 LET LO#eL@-F*x9 

230 NEXT X 

240 FOR Y2Y¥!+vY0/2 TO Ye-¥Y8/2 STEP YB 
258 LET FeeaKOal axl /(X 14K 1+Y*Y) 
26@ LET LO=LO-F*y¥@Q 

270 LET FxeaK Gal aX2/(X 2eX O4+Y*Y) 
286 LET LO2zL@+FayYQ 

290 NEXT Y 

3008 PRINT “LINE INTEGRAL #"3L@ 
318 PRINT 

326 GOTO 168 

338 END 


Since F = (F1,F2,F3) is due to a long straight wire (uol/27r, tangential 
with sense by the right-hand rule), you expect §Fodk to equal zero unless 
the rectangular contour encloses the z-axis. When the contour encloses the 
Z-axis, you expect $Fed& = wol in mks units. A run of this line integral 
program looks like: 


RUN 
AMPFAR 


(XsY) FOR 2 CORNERS? 

Pe lsneds lel 

LINE INTEGRAL 2 1-25664E-@6 
(XsY) FOR 2 CORNERS? 

7 la ls 252 

LINE INTEGRAL 2=-3-77476E-15 


(XeY) FOR 2@ CORNERS? 
? 


Practice Exercise: Modify the program for two, long straight wires parallel 
to the z-axis, carrying currents of +7 and -1, and cutting the x-y plane at 
(+.5,0) and (-.5,0), respectively. Where should the line integral be zero and 
where non-zero? 


Practice Exercise: Modify the program for a single line charge coinciding 
with the z-axis and having a linear charge density, \, of +1. Where should 
the line integral be zero and where non-zero? 


Further examples are considered in the problems at the end of the chapter. 
—_> 
The vector field F could be given as equations or as data. 


CONCLUSION 


This chapter has displayed pictures of several vector fields. You were asked 
to try to identify which pictures could represent electrostatic and mag- 
netostatic fields. The chapter then showed you a way to compute the 
surface and line integrals necessary for Gauss’s, Ampere’s and Faraday’s 
Laws in order to check vector fields quantitatively. By working a few 
examples from the exercises, you can deepen your understanding further. 





CHAPTER FOUR EXERCISES 


1. Consider the field distribution created by four line charges parallel to 


the z-axis. The line charge densities are: +1 at (1,1), +1 at (1, -1), -1 at 
(-1,1) and —1 at (-1,-1). 


a} Consider the region around the (1,1) line. What is the surface 
integral? Is it dependent on the size of the box around the charge? If 
so, is it dependent on all lengths or just one? What is the line integral 
around this point ($F +d)? Is it dependent on the size of the 
contour? 


b}) Consider the region of space for positive x that includes both the 
lines at (1,1) and (1,~1). What is the surface integral ($F dA)? Is it 
dependent on the size of the _box? If so, on what does it depend? 
What is the line integral (¢F° di) around this region? Is it dependent 
on the size of the contour? 


~aa 


c) Consider the region including all the line charges. What is the surface 
integral ($F *dA)? Is it dependent on the size of the box? If so, how 
is it dependent? What is the line integral (¢F°dX) around the region? 


ls it dependent on the size of the contour? 


d) Consider the region inside the four points. What is the surface 
integral (¢F*dA)? What is the line integral ($F *d&)? 


e) Do your results agree with Gauss’s and Faraday’s Laws? 


. Consider the magnetostatic field distribution created by four current 
lines parallel to the z-axis. The positions and currents are: +1 at (1,1), 
+1 at (1,-1), -1 at (-1,1) and -1 at (-1,-1). 


a) Consider the region around the (1,1) line. What is the surface 
integral (¢F*dA)? Is it dependent on the size of the box around the 
charge? If so, is it dependent on all lengths or just one? What is the 
line integral around this point (¢F +dXJ)? Is it dependent on the size of 
the contour? 


b) Consider the region of space for positive x that includes both the 
lines at (1,1) and (1,-1). What is the surface integral ($F dA)? Is it 
dependent on the size of the box? If so, on what does it depend? 
What is the line integral ($F d2) around this region? Is it dependent 
on the size of the contour? 


ee 


Consider the region including all the line charges. What is the surface 
integral ($F dA}? Is it dependent on the size of the box? If so, how 
is it dependent? What is the line integral ($F de) around the region? 
ls it dependent on the size of the contour? 


Cc 


d) Consider the _region inside the four points. What is the surface 
integral ($F dA)? What is the line integral ($F AR)? 


e} Do your results agree with Gauss’s and Faraday’s Laws? 
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3. Consider a combined electrostatic and magnetostatic field distribution. 
Part of it is created by two line charges parallel to the z-axis: +3 linear 
charge density at (1,1) and —3 linear charge density at (-1,1). The other 
part is created by two fine currents parallel to the z-axis with currents 
and positions: +2 at (1,-1) and —2 at (-1,-1). 


a) Consider the _fegion around the (7,1) line. What is the surface 
integral (§F*dA)? Is it dependent on the size of the box around the 
charge? If so, is it dependent on all lengths or just one? What is the 
line integral around this point ($F +dQ)? Is it dependent on the size of 
the contour? 


b) Consider the region of space for positive x that includes both the 
lines at (1,1) and (1,-1). What is the surface integral ($F +*dA)? Is it 
dependent on the size of the box? If so, on what does it depend? 
What is the line integral (¢F+d2) around this region? Is it dependent 
on the size of the contour? 


c) Consider the region including all the line charges. What is the surface 
integral ($F *dA)? Is it dependent on the size of the box? If so, how 
is it dependent? What is the line integral {$F «d) around the region? 


ls it dependent on the size of the contour? 


d) Consider the region inside the four points. What is the surface 
integral (6F*dA)? What is the line integral ($F +9)? 


e) Do your results agree with Gauss’s and Faraday’s Laws? 


4. Consider the following field. 


P= 0 
FY =0 
F,=1forz>0 
F,=0Oforz=0 
F,=-1forz<0 


a) What is the surface integral for the box defined by (-1,- 1,-1) and 
(1,1,1)? lf it is not zero, is the surface integral dependent on the size 
of the box centered at (0,0,0)? In what way? 


b) If the surface integral is not zero, can you identify the charge 
distribution that creates this field? 


c) What are the line integrals for the square loops defined by (1,1,0) 
and (-—1,-1,0); (1,0,1) and (-1,0,-1), (0,1,1) and (0,-1,-1). 


d) If the line integral is not zero, can you identify the current distri- 
bution that creates this field? 
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5. Consider the following field. 


f= 0 
FY =0 
F,=1forx >0 
F,=Oforx=0 


F,=-1forx<0 


a) What is the surface integral for the box defined by (-—1,-1,-1) and 
(1,1,1}? If it is not zero, is the surface integral dependent on the size 
of the box centered at (0,0,0)? In what way? 


b) If the surface integral is not zero, can you identify the charge 
distribution that creates this field? 


c) What are the line integrals for the square loops defined by (1,1,0) 
and (~1,~1,0),; (1,0,1) and (-1,0,-1); (0,1,1) and (0,-1,-1). 


d) If the line integral is not zero, can you identify the current distri- 
bution that creates this field? 


6. Consider the following field: 
> A 
Fer * r(2tr?) 
a) What is the surface integral for the box defined by (-—1,-1,-1) and 
(1,1,1)? If it is non-zero, is the surface integral dependent on the size 


of the box centered at (0,0,0)? In what way? 


b) If the surface integral is not zero, can you identify the charge 
distribution that creates this field? 


c) Can you prove that the line integral for any field that can be written 
“Aw 
as F = rf(r) is zero? 


7. Not ail fields are electric or magnetic. For example, it is often worth- 
while to consider velocity fields of fluids. Consider the following fluid 


velocity field: 
V,=0 
Vy =Q 
Vy =(1-y7/4) 


a) What is the surface integral for the box defined by (-1,-1,-1} and 
(1,1,1)? lf it is non-zero, is the surface integral dependent on the size 
of the box centered at (0,0,0)? In what way? 


b) What does this surface integral mean physically {i.e., are there 
sources or sinks for the fluid within the region)? 
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c) What are the line integrals for the square loops defined by (1,1,0) 
and (-1,-1,0),; (1,0,1) and (-1,0,-1); (0,1,1) and (0,-1,-1)? 


d) What is the line integral around the square loop defined by (0,0,0) 
and (1,1,0)? 


CHAPTER FIVE: THE LAPLACE AND 
POISSON EQUATIONS 


INTRODUCTION 


A combination of Laplace’s and Poisson’s equations represents one of the 
fundamental ways to describe electrostatic potentials. (Similar equations 
arise in many other branches of e & m and of physics in general.) You 
sometimes know the potential, V, everywhere on some surface that encloses 
a definite volume, and you want to find the potential throughout that vol- 
ume. If the volume has a known volume charge density, o(x,y,z) (measured 
as charge per unit volume), then you can calculate the potential everywhere 
using Poisson’s equation. If the volume charge density, p is zero through 
the region, then Poisson’s equation reduces to Laplace’s equation, one of 
the simplest partial differential equations used commonly in physics. 


Poisson’s equation in three dimensions is 


Vv eV, dV 
Me ee ee Amko p(x,y,2) (15) 





When p(x,y,z) = 0 you have Laplace’s equation 


07 87 O0?V_ 
xe tage + 33 =0 (16) 


Derivations of these equations and discussions of their analytical solutions 
can be found in standard physics textbooks. 


These equations can be solved easily on a computer, and they represent 
simple examples of how partial differential equations can be solved nu- 
merically. For convenience in displaying results, we will consider two 
dimensional cases involving only x and y. The method of solution is valid 
for three dimensions, too, but more complicated partial differential 
equations sometimes demand more sophisticated methods of numerical 
solution than the method we will present. 


DIFFERENCE EQUATIONS 


The simplest method for solving partial differential equations numerically 
is based on replacing the differential equation (involving infinitesimal dx, 
dy, dz) by a difference equation (involving finite differences Ax, Ay, Az). 
There are three criteria that the resulting difference equation must satisfy: 
1) it must represent the partial differential equation correctly in the limit 
as Ax>0, Ay>0O and Az~>0; 2) it must be stable (which means that, if a 
small error occurs at some numerical step, the error must not grow); and 
3) it must converge to some value (in other words, it must produce an 
answer). Laplace’s equation is such that a simple difference equation 
approach can be shown to satisfy al! three requirements. If the volume 
charge density, 0, does not vary too wildly in space, the method will also 
work for Poisson’s equation. 
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We start by defining a grid of points, labelled by i and j (for our two 
dimensiona! problems), covering the region where we want to compute the 
potential, V. How would the first partial derivative, 0V/dx, be written in 
terms of the values of Vij on the grid points? You could use the forward 
difference expression for a derivative 0V/dx~(Vj44_ j — Vjj)/Ax or you 
could use the corresponding backward derivative. A better approximation 
is called the central difference approximation and uses grid points in front 
of and behind the grid point (i,j) in question. Then 0V/dx*(Vj+1,j — 
Vj-4 j)/(2Ax). 


The second partial derivative, 0*V/0x?, equals 0(0V/dx)/dx so we can 
apply the same idea of approximating the partial derivatives by difference 
quotients. Using the central difference approximation, 0? V/dx? =(Vi44 = 
2Vii ij + Vi-4 pdx)? . 


Practice Exercise: Show that the difference equation approximating 
07 V/ax* comes from two applications of the central difference approxi- 
mation for the first derivative. The first application gives the first deri- 
vatives at (i+.5,j) and (i-.5,j); the second application yields the second 
derivative at (i,j). 


THE DIFFERENCE EQUATIONS FOR LAPLACE’S AND 
POISSON’S EQUATIONS 


Using the expression for approximations to the second partial derivatives 
of V on the grid of points, the (two-dimensional) Laplace’s equation 
becomes 


Vier + Vina + Viger + Vij-1 - AV = 9 7) 
assuming Ax = Ay. 
Practice Exercise: Derive this difference equation for Laplace’s equation. 


Practice Exercise: Derive the three-dimensional difference equation for 
Laplace’s equation. 


This difference equation can be solved for the potential at the grid point 


(i,j): 
Vij = (Vi4r,j + Vi-1j + Vij+1 * Vi,j-1)/4 | (18) 


which says the potential at (i,j) is the average of the potentials at the nearest 
neighbor grid points. This equation for V suggests a simple strategy on the 
computer: Start with any values of the potential inside the region and with 
the known values on the grid points on the surface; moving sequentially 
through the interior grid points, replace the inside values by the average of 
nearest neighbor values; then repeat the process until interior values of the 
potential stop changing. The boundary values, which remain constant, will 
gradually influence the values throughout the interior. When all the inside 
values become equal to the average of their nearest neighbors (that is,after 





all the values of the potential at interior points stop changing), then you 
have found the solution of Laplace’s equation with the given potentials on 
the surface. 


Poisson’s equation is just as simple. The difference equation (in two 
dimensions) becomes 


(Vier jt Viet + Vijer + Vij-1 7 4Vi,j) = Rij (AS) (19) 


where Ax = Ay = AS and where Rij is just the (known) quantity, 
Ark op(x,y), at the grid point (i,j). Ko sets the units. Since Rj | does not 
change value at each particular (i,j) as you iterate over the grid, programs 
solving Poisson’s equation by the difference equation approach look very 


similar to those solving Laplace’s equation. 


Example 1. Find the potential V in a charge-free, two-dimensional, square 
region given that the potential varies linearly around the edge from V =0 
at one corner to V = 18 at the opposite corner. 





A program implementing the iterative strategy to solve Laplace’s equation 
for this problem might look like the following: 


LAPLAC 


160 DIM VC 18s161,U0 18. 18) 

11@ READ VC 121],V0 18.1]. VC 1G, 181], VE 1s 18) 
128 DATA 029,18.9 

136 LET VOeVC ltl, 13+VC 1, 1034+ VE 18. 18)+V0 18, 1) 
14@ FOR I=2 TO 9 

156 LET V€1l,11]=I-1 

160 LET. V016,1129+1<1 


17@ LET VCls1]8I-1 Abesied 
18@ LET VE1Is161=9+I-1 Hapa 
192 LET VOsVO+VE lo LI+VCIG, 13+ VCI. 134+VC1, 18) 

200 NEXT 1 

210 FOR 122 TO 9 

228 FOR J#2@ TO 9 tnitialize interior Values 

230 LET ViI».3]=VO/36| to Average of 

240 NEXT J Boundary Values 


27@ FOR J=1 TO 18 

28@ LET ULLsJ}=VElsd} | pe is 
298 NEXT J 

300 NEXT I 

31@ LET Vise 

326 FOR 122 TO 9 


338 FOR J=2 TO 9 fei 
340 LET Vi lsd]=CUCI + bod] FULT = Led I4UC Lo d+ LI +UL Lod-13) 74] porenrials b, 
35G IF ABSCCVEL»JI-UCIsJ}I>/VEIsJ3)<Vl THEN 372 Averaging 
368 LET VISABS(C VEL, JI-UCIsd))/VET2II) |Ynew > Vora Neighbor 
370 NEXT J V Values 


new 


38@ NEXT I 

398 LET M=M+1 

AGS IF Vl>-@@1 THEN 260 -— Jest for.1% Accuracy 

416 PRINT "# OF ITERATIONS ="3M 

42@ FOR J=1@ TO 1 STEP -1 Print bie 
430 FOR I=] TO 164 Potential 
446 PRINT Vi IlsJI)5 Matrix 

45@ NEXT I 

4686 PRINT 

470 NEXT J 

486 END 


The program initializes the potential at interior points of the region to the 
average of the values on the boundary. It is a property of Laplace's 
equation that the absolute maximum and absolute minimum values of 
potential must lie on the boundary. Initializing the interior points to the 
average value saves some time in the calculation. 
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It turns out that Laplace’s equation has the interesting (and unusual) 
property that, when you average the nearest neighbor potentials, you can 
use some old values (from the last iteration over the grid) and some new 
values (from the present iteration). This allows you to use only one storage 
matrix for the potential. | 


Practice Exercise: Modify the program so that you use only the matrix V 
and average two old and two new nearest neighbor potentials, 


A RUN of the program given above looks like: 





RUN 
LAPLAC 
# OF ITERATIONS = 34 
9 16 11 12 13 
14 15 16 17 18 
8 Se 9.99882 180996 11.992 
12.9885 13.9868 14.9885 1529932 17 
7 Be BB119 9.98801 9.99541 19.9885 
1169815 12.9783 13.9893 14.9885 16 
6 7 «88402 $-G8461 9-60601 9299137 
18.9823 11-9767 12.9783 13-9868 15 
5 6¢ 68882 7681156 8.908 65 9.88001 
9.98955 10.9823 11-9815 12.9885 14 
4 5.81156 6.81848 7681771 8°@1047 
9.80001 9099137 10-9885 11-992 13 
3 4.81325 5.82172 6.8233 7.81771 
8-808 65 9.99681 9.99541 180996 12 
2 3681155 4e81967 5.82172 6281848 
7081156 Se OQ46i 9.88901 9.99882 13 
1 2.88683 3-81155 4.81325 5.81156 
6-298 G2 798401 8.88119 96 12 
g ! 2 3 4 
5 6 7 8 9 


Example 2. Find the potential V throughout a square region, given that 
the potential on the boundary varies linearly from 0 at one corner to 18 at 
the opposite corner and that a uniform surface charge density of +.5€ 
coulombs per m? fills the region. 





The problem calls for the solution of Poisson’s equation in a situation very 
similar to that above for Laplace’s equation. A program that solves this 
Poisson equation situation might look like: 


POISSO 


160 DIM V0C16.181,UE 18.18) 

118 READ VC ls.1).V0 18. 1)5 V0 18s 101]5 VC 1, 18) 
128 DATA 0.9. 18,9 

138 LET V@=VC ls. 1)4+VC 1s 1614+VE 10s 1803+VC 18, 1] 
149 FOR I=2 TO 9 

1S8@ LET VOl,zIl]#1-1 

168 LET VC 1l1@s1) 29¢!I-1 

178 LET VCI,-l}=1I-) 

188 LET VUIs 181=9+Il-1 

190 LET VOeVO+VEC ls I1}4VC1G,134+V0I, 134+VC 1s 183 
260 NEXT ! 

210 FOR I=2 TO 9 

220 FOR J=2 TO 9 

230 LET Vf1,J12# V8/36 

246 NEXT J 

28@ NEXT I 

268 FOR I#1 TO 180 

270 FOR J=!1 TO 18 

26e@ LET UCIsJI=VCI.JI] 

296 NEXT J 


380 NEXT I 

318 LET Vl=@ 

328 FOR I#2 TO 9 

3380 FOR J#2@ TO 9 

348 LET R-5 

356 LET VOI“JIACUCI + ls JIFUCI <1, 5I+UCI,J+13+UC Ls J-13+R) 74 
360 IF ABSC (CVE I,JI“ULIZJID/VLIZ,JI)<V1I THEN 388 
376 LET VISABSCCVCIs JI] -UCILJI) /VEI.II1) 

388 NEXT J 

398 NEXT I 

406 LET M=M+!1 

4109 IF V1>.881 THEN 268 

420 PRINT “*# OF ITERATIONS ="3M 

43@ FOR J21@ TO 1 STEP ~1 

446 FOR I=] TO 18 

458 PRINT VCIsJI]> 

468 NEXT I 


478 PRINT 
482 NEXT J 
49@ END 


Practice Exercise: Annotate the Poisson program by written comments by 
each line. 


A run of this program looks like: 


RUN 
POISSO 
# OF ITERATIONS = 48 
16 11 le 13 
14 1S 16 17 18 
8 9658738 18927 12-1182 13-2948 
14-6 2044 15-1172 15.9259 16-58 66 17 
7 8292723 1@e5112 11-8531 13-8106 
14.8899 14-8513 15-589 1 1509259 16 
6 Be 11986 985412 11.2951 12¢5807 
13-4998 14-2926 14-8513 15-1172 i5 
5 7220652 99-8131 18+«5626 11.7922 
1207311 13-4998 14-8899 1422844 14 
4 6- 2869 8+81384 9+58357 18°7333 
11.7322 12-5687 13+¢8106 13-2048 13 
3 Se 12085 6685596 8.29763 950357 
18-5626 11-2951 11-8531 12-1182 12 
2 3-9 2835 §-51334 6685596 8-@1384 
9-@813i 9e85412 186-5112 18-927 11 
1 2-58811 3¢9 2835 5-12825 6-20 69 
720652 Ge 11986 892723 9-58 738 12 
@ 1 2 3 4 
5 6 7 8 9 


Notice that the final interior values have changed from those of the 
Laplace solution. The symmetries have been maintained because the 
charge density is uniform. 


Practice Exercise: Identify lines of symmetry in the values of potential. 


Terminal plotting can also be used in a manner similar to the terminal 
plotting of potential done in Chapter One. 


Practice Exercise: Modify the Laplace and Poisson programs to print out 
characters representing the values of the potential at each point on the 
grid. 
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CONCLUSION 
The Laplace and Poisson equations are very basic to electrostatics. This 
chapter showed you one way to solve these equations with the computer. 


The numerical method of solution presented is also applicable to other 
partial differential equations arising in physics. 


CHAPTER FIVE EXERCISES 


1. Consider the square region with the boundary potentials: 


0 2 0 
2 4 
4 8 
6 12 
8 16 
10 20 
8 16 
6 12 
4 8 
2 4 
0 5 0 


(constant zero potential on two sides and increasing linearly to the 
midpoint then decreasing to zero again on the other sides). 


a) Determine the potential inside the region. 


b) From a character plot of the potential in the region, draw several 
equipotentials. 


c) Describe the symmetries in the potential. 


2. With the boundary potential described in Exercise 1, calculate the 
potential in the region when there is a uniform charge density of —.5€5 
coulombs/m? throughout the region. 


a) From a character plot of the potential in the region, draw several 
equipotentials, 


b) Describe the symmetries in the potential. 


3. Consider the following line charges parallel to the z-axis: +1 at (6,0); +1 
at (-6,0); +1 at (0,6); and +1 at (0,-6). 


a} Using the method discussed in Chapter 1, determine the potential in 
the square defined by the corners [5,5] and [-—5,-5]. (Find the 
value of the potential on a grid with ten points on each side.) 


b) Use the values of the potential obtained in Part (a) for the potentials 
on the edges of the square region defined by the corners [5,5] and 
[-5,-5]. Use Laplace’s equation and determine the potential inside 
this region. 





c) Compare the results obtained by the two methods (Parts (a) and 
(b)) for interior points. 


. The symmetries observed in the potential pattern in the text and 
Exercise 1 can be used to cut down the number of points at which the 
potential needs to be evaluated. For any symmetric pattern, we only 
need to calculate the potential for a smaller region of space and then 
use the symmetry to copy the potential into the other regions. 


Consider the potential in Exercise 1. If you consider the pattern 
centered at the origin the pattern above the x-axis is mirrored below the 
x-axis. To solve the problem using the symmetry, set the potential the 
same way for positive y but when you solve for the potential along the 
x-axis (the line of symmetry), consider the value of each point just 
below the axis to be the same as the value of the point just above the 
axis. The grid point equation for this line of symmetry becomes: 


Vij = (Vii4 i + Vi_4 sj + 2V js j-1)/4. 


a) Resolve the potential for Exercise 1 taking advantage of the 
symmetry. 


b) Resolve the potential in Example 2 in the text using the symmetry 
along the diagonal, 


. Consider the square region with the potential along the edges defined 
by: 
24 68108 6 4 2 


0 0 

2 2 
4 4 

6 6 

8 8 
10 10 
8 8 
6 6 

4 4 
2 2 

0 0 


246 8108 6 4 2 


(with all edges identical, with the corners at zero potentials, and the 
potential increasing linearly in from the corners to the midpoint where 
the potential is ten.) 


a) What are the symmetries? 


b) Taking advantage of all of the symmetries, calculate the potential 
inside the region. Make a character plot of the potential values. 


c) Draw in some of the equipotentials, 


. Consider the same square region as Exercise 4 but with a uniform 
charge density of .5€ (|x| + ly|) coulombs per m? throughout the 
region. 
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a) What are the symmetries? 


b) Taking advantage of all of the symmetries, calculate the potential 
inside the region. 


c) Make a character plot of the full region and draw the equipotential 
lines. 


. Consider the following potential: a square boundary with a potential of 
zero on it. Let the square be defined by the corners (-9,-9) and (9,9). 
Inside the boundary is a square defined by the corners (-—3,-3) and 
(3,3) at a constant potential of 100. The potential in the region in 
between can be solved by the techniques described in the text. You 
must check to see if the point being evaluated is one of the center 
points with a fixed potential and, if it is, go on to the next point leaving 
that grid point’s potential at 100 as initially defined. 





a) What are the symmetries in this problem? 


b) Using the symmetries, determine the value of the potential inside the 
region and make a character plot of the full region. 


c}) Draw in some of the equipotentials. 


. Consider the potential caused by a +1 and —1 point charge enclosed in a 
square whose edges are held at zero potential. In order to approximate 
this situation, assume that the potential near the charge is constant and 
is unaffected by the rest of the region (i.e., assume the potential near 
the +1 charge is determined by the equation kgq/r). 





Let the corners of the zero potential box be defined by (-5,-—5) and 
(5,5) and the plus and minus charge positions: +1 at (+2.5,0) and —1 at 
(-—2.5,0). 





a) What are the symmetries in the problem? 
b) What can you say about the potential on the y-axis? 


c) How can the solution of Part (b) be used to simplify this problem? 


d) Compute the potential and make an expanded character plot of the- 


potential. 
e) Draw in the equipotential lines. 


f) Compare your results in shape to the simple dipole charge 
distribution. 


LAPLACE AND POISSON EQUATIONS () 67 








68 O ELECTRIC AND MAGNETIC FIELDS 


APPENDIX A: TERMINAL PLOTTING 


The subroutines in TTYPLO (listed below) allow you to plot a set of 
curves on your terminal. The subroutines use the letter O variables, so you 
should not use these variables in your programs. The output from the 
subroutines is a 5’ by 5” x-y plot (assuming 6 lines/inch and 10 
characters/inch) with the x-axis across the page and the y-axis up the page. 


The range of x values on the plot are set by defining values for X8 = the 
left side of the plot, X9 = the right side of the plot, Y8 = the bottom of 
the plot, and Y9 = the top of the plot. The subroutines are called with 
“GOSUB 9000”. 


To plot a point on the graph, you define the values of X@ = the x coordi- 
nate and Y@ = the y coordinate of the point. Also, specify a value for ZQ, 
which will be the character plotted at the point (X@,Y@) on the graph. 
You then call ‘‘GOSUB 9190”. 


The allowed values of Z@ are O through 9. Since these variables 
X8,X9,Y8,Y9,X0,Y0,Z@ are used in these special ways for the plotting, 
you should be careful how you use them in your programs. 


To produce the final plot on your terminal, call “GOSUB 9200’’. A sample 
of the use of TTYPLO follows. The program plots a sine curve (as the 
character “’1’’). 


SINE 


180 READ X8sX9sY8sY9sZ98 

118 DATA Gs 6-28318s.-14 151 

128 GOSUB 9282 

13@ FOR X@=X8 TO X9 STEP (X9-X8) /58 
146 LET Y@=SIN(CX@) 

1586 GOSUB 9186 

160 NEXT X6 

176 GOSUB 9288 

188 STOP 
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APP- TTY PLO 
RUN 
SINE 


MAX Y= 1 
XXXXXXXXXXKXKXXXXXXXXXKKAKXKKKKAKKXK AXRKXKAXKXXAXKXKXXXAXX XX XX 
rlil 
1} il 


Sem KKM KR K RR KKK 


b111 
XXXXXXXXXAXKKXXAXKXKKXAKAKXAKARKKKXAAKAAKAKKX AK KKXAKKAKKAKKK 
MIN Ys~1 MIN X= 6 MAX X= 6628318 


END 


90880 REM INITIALIZE PLOT 

9819 DIM 0£366] 

9220 FOR O01i=! TO 306 

9030 LET OCO11]=¢ 

9042 NEXT O01 

9@58 LET 0(523=252 

9260 RETURN 

91@@ REM STORE POINTS 

9105 IF (X@-X8)*(X@G-X9)>6 THEN 9165 
9110 IF CYO-Y8)*(Y@-Y9)>@ THEN 9165 
9115S IF Z@ >= @ THEN 9125 

9126 LET 26=¢ 

9125 IF ZO <= 9 THEN 9135 

9130 LET Z2@=9 

9135 LET O0€52)=0(652]+1 

9149 LET OC OC523]=100Ge INTC 38% (YO-Y8) /(Y9-Y8)+1-5) 
9145 LET Of 0£5271=0C€ OF 5212+ 1G" INTC SO*(XG-XB) /(X9-XB + 1059428 
9150 IF 0€52)<396 THEN 9165 

9155 PRINT “PART OF THE PICTURE NOW" 
916@ GOTO 9228¢ 

9165 RETURN 

9200 REM PRINT OUT PLOT 

9218 PRINT "MAX Y="3_Y9 

92289 PRINT “ XXXXXXXXXXXXXXXXXXXXXXXXXKX KKK XAXKKX ANAK KK AKKKK KXKK KX 
9236 FOR 0283! TO 1 STEP -1 

9248 FOR 0321 TO 51 

9258 LET 0£03}=-1 

9260 NEXT 03 

9278 LET 06=¢ 

928@ FOR 03253 TO 0C52) 

9298 IF 0O£031<@ THEN 9362 

9380 IF INTC0£03]/19800) <> O2 THEN 9366 
9319 LET O1=0€03)-1800*«INTC Of 032 /1090) 
9320 LET OCINTCO1/1@)2=#01-10"INT( 61/18) 
9330 IF O6>INTCO1/1@) THEN 9350 

9348 LET O6"INTC(O1/19) 

9358 LET 0(03)=-1 

9368 NEXT 03 

9370 PRINT "Y"3 

9380 FOR 03=1 TO 06 





9398 GOTO 0£F033+2 OF 94390,9428,9 446.9 468.948 6,95002952859540,9560s9 588. 


9598 


9486 
9418 
9428 
9438 
9449 
9458 
9468 
9478 
9482 
9496 
9588 
9516 
9526 
9538 
9548 
9558 
9568 
9576 
9588 
9598 
96088 
9616 
9628 
9638 
9648 
9658 
9668 
9676 
9689 


PRINT * ‘3 
GOTO 9616 
PRINT "@°*3 
GOTO 9618 
PRINT "13 
GOTO 9618 
PRINT “2's 
GOTO 9618 
PRINT "3°3 
GOTO 9618 
PRINT "&4°*3 
GOTO 9618 
PRINT "53 
GOTO 9616 
PRINT "6" 
GOTO 9616 
PRINT “73 
GOTO 9618 
PRINT "6 °3 
GOTO 9618 
PRINT "9" 
NEXT 03 
PRINT 
NEXT 02 
PRINT ™ 


PRINT "MIN Y="3YS3 °° 


GOSUB 9888 
RETURN 
END 


XXXXAXAXAMRKK AKAMA A AA X ** 
MIN X="3X6, MAX X="3 X9 
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APPENDIX B: THE HALF-STEP 
INTERATIVE INTEGRATION 


Numerical methods comprise an entire subject in mathematics. For the 
moment, all we need to show you is that the half-step method used heavily 
in the text would seem to give better answers than the simplest approxi- 
mation (which is called Euler’s method). The problem we are discussing is 
a general one. A physical law often can be stated in terms of the derivative 
of a function you wish to find; you integrate the differential equation to 
find the answer. 


Consider a general curve for y=f(x); suppose you know the value of y at 
X=Xq and you have a way to calculate the derivative of y with respect to x 
anywhere. The problem is to get as good an approximation to the correct 
value y=f(x9+Ax) as possible. The simplest method (Euler's method) uses 
the derivative at Xo and approximates f(x g+Ax) as 
£(x9 tAx)=F(x9 )+#’(xq JAx. 


Euler’s method is shown on the sketch by the dashed line. Clearly, the 
method is correct in the limit as Ax-?0. For finite Ax the method 
produces excessive errors. 


(x) 
+ 


Euler approximation 
a 


“a 


true 
f(x) 





x x + AX x—> 


The geometrical interpretation of the theorem mathematicians call the 
Mean Value Theorem says that there is some point on the (continuously 
differentiable) curve y=f(x) such that the derivative at that point has the 
same value as the slope of the chord from (xo,f(x9)) to (x9 tAx,f(x9 +Ax)). 
If the theorem told us where that point was, everything would be easy. 
So-called higher order methods of iterative integration (such as the fourth- 
order Runge-Kutta method or predictor-corrector method) are better and 
better ways to approximate this ‘‘correct’’ value of slope. 


Our half-step method is based on the fact that, in general, evaluating the 
derivative near the center of the interval [xo,X9t+tAx] is better than using a 
derivative at the end of the interval. In fact, you can show that, if f(x) is a 
parabola, then a true half-step method is exact. This means that, if you 
expand f(x) around xg in a Taylor series, the half-step method will be 
correct through terms in (Ax). That is why half-step methods are called 
“‘second-order’’ Runge-Kutta methods. 
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In several cases used in the text of this unit an approximate half-step 
method, one that approximates x9+Ax/2 by using the Ax calculated at the 
last step, is used instead of a true half-step. Practically always, such an 
approximate xj+Ax/2 lies near the center of the interval [xg,x9t+tAx] and 
actually increases the accuracy of the method. 


Practice Exercise: Using any of the iterative integration programs from the 
text, compare Euler’s method (which uses the derivative at the point xq) 
to the half-step method. Show that, for any given step size, the half-step 
method is more accurate and that, to achieve the same accuracy as the 
half-step method, a much smaller step size (and hence many more iter- 
ations) must be used by Euler’s method. 
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ANSWERS TO SELECTED 


EXERCISES 


CHAPTER ONE 


Exercise 1. 


KeY>) OF LERT-BOTTON & RIGHT= TOP? 


T =6S52—-205, 7 eo 4$e 205 

TNs & MAX, VT 

? -100. <1! 

SZ ZY VV YX XXX AXAKAXKXKKKAAXV VV VY S22 
AVY VYXXXX SUV OWE XXXXYYYYZ 
YVVEXX WS VUV VY SUS VV VUE WUXXXT TTY 
YURXXWEVVUSUT TTT TT TUVUUVY UWXXXTY 
TAX WS VVEUT SSRAGAORRS STUUVY UUXXY 
RAVWVVUTSROPOGUNNGS PERS TY VY UVXX 
MVVVUTSRPOMK HHH I KMO PRS TUVY UXX 
WEVUTSRENLI EABSBALI LY PRSTUVYUUX 
WVVUTSQGN! BE 0000GE DING OSTUYV UWX 
WMVVUTROOLH B! 00000 1 BHLOGRTUVY WX 
WVVVUTS OOM! DE0 00006 DING ESTUVVWX 
WUVEUTSRPVULI EARSTARILN PRSTUUW UX 
RKVVV UT SRPOMK JHHHIKNG PRS TUVV UXX 
XKVEVUUT SREOPOONNNE SPARS TUVV WUXX 
VAX VV VU UUTSS RREGERASSTUYUVVEVXXY 
YUAX AWW VURUE TT TTT? TUG UY BUXXK YY 
LV VRAX PU VS VVC UY VY UUUXXXYYTY 
RVVVYXXXX WU VOW GUN EW NX XXKYY YYZ 
SEZ VV VV VAR XXX AXAKAAXXAAKVVVVY ZZE 


KoY) OF LEPT+BOTTON & RIGHT-TOP? 


? 045, ¢4 0S. « $$, « O$ 

MMe &2 WAX. V? 

? 402200 

O1 1129233334448444443393392221110 
1212223344455 5666663555 443322211 
1223344556999 SESE6 PRP E665 443322 1 
223445667899 AABBBAAP 9S 266544322 
23 445676 9 ABS DEEP EE DE BASE 76535433 
3445679 ABDFE! JKLKS 3 GF DBAS 365543 
345679 AG EBM PSUUUS PN J 6 E6AD 76543 
4667S AGESK OTEZLZZELL TOXS ECA 3654 
GETS ADFINSEZLZLZZZLSNI FDASTES$4 
46679 BOG OUZ 222222200) 6 DB? 3654 
46679 ADPEMSEZZLIZIZ LE SN F DAD 7654 
@@OFE AGEGKOTZELZZZEZTOKGEC AS 76S4 
345679 ASEGIH PSUUYS Ft 6 ESAS 76543 
3443679 ABDF OI UKLKJ 1 65 BBAD 765543 
234486769 ABCDEL? EE DEBADE 76539433 
223445667699 AABBBAAS 96 ¥66544332 
12233445 3637768888 F77665443322 1 
112223344455 $666665555 443322211 
Ol § 12223333 4444444443333 2221110 


CKXeoY) OF LEFT+BOTTOM & RIGHT<TOPT 


? -$,-$8,5.5 

MEN. & MAXe VP 

? 6.5 
9222292222233333333333333322222 
2222222233333333333333333333322 
2232222393933333344444443333933332 
922223339333344044444444443933333 
22223333334444555558555444843333 
2228333334448 8666986665 39444333 
22239333344458968699968E 7665844433 
282333334445679 BEEDGAS7768 54443 
92339333344444494 TUOH BADF6S $4443 
#2333333 4443000SE221.FB9 76684443 
92233333 344444494 TVGH DAS 365 54443 
2823333344456 79 BEEDCAS F76554443 
222333334448 867899988 7665544433 
223239333 4445566677666558444333 
22223333933444455588555544643333 
22292333333344464444444444333333 
9822222333333334444844333333332 
222292223333333393333339339333322 
22 223222222333333933393333322222 
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f) Near enough to each line charge, all other charges can be neglected 
and the equipotential is that of a single line charge (circular). Far 
away from all three tine charges, the equipotential lines are again 
circular and appear in value as if they were created by a single line 
charge with a linear charge density of +3. 


Exercise 3. 





f) The x-z and y-z planes are the zero potential planes. 


Exercise 4. 


c) Near enough to each line charge, all other charges can be neglected 
and the equipotential is that of a single line charge (circular). Far 
away from all three tine charges, the equipotential lines are again 
circular and appear in value as if they were created by a single line 
charge with a linear charge density of +4. 
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Exercise 5. 


x 


LX 





c) The point charge model of a tine charge and the true line charge 
differ in three places: (1) at the ends of the nine charges, (2) very 
close to the line, and (3) very far away from the line segment. 


d) To make a better approximation, increase the number of charge 
points on the line and make the point charge line segment longer. 


Exercise 6. 


c) The plane defined by alli points with x = .5 and the plane defined by 
all points with x = -.5 have zero potential. 


Exercise 7. 


a) A zero potential plane on the y-z plane and the charges +1 at (.5,.5) 
and -1 at (~.5,+.5). 


b) The x-z and y-z planes. The charge is +1 at (.5,.5). 





c) Exercise 6 is (an approximation to) the potential for a point charge 
at the origin and two conducting planes of zero potential defined by 
the set of prints with x = ~.5 and the set of points with x = .5. 


d) Solve this by using the following charge distribution: The original 
charges and the image charges -2 at (-.5,0); +1 at (~-.5,.5) and +1 
(-.5,-.5). 


Exercise 8. 
c) The potential inside the cylinder is constant. 


d) Gauss’s Law is $E°dA = 4nkgq. Since any surface constructed inside 
the surface will contain no charge, the field is zero, Outside the 
cylinder, the field is the same as that of a single line charge with a 
linear charge density equal to the sum of the line charges used to 
make the cylinder. 


CHAPTER TWO 
Exercise 1. 


(3, 3) 


(~3, -3) 


d) Close to each wire, the field is similar to a single current wire of that 
wire’s current. Far away, the field is similar to a single current- 
carrying wire with a current given by the (algebraic) sum of the 
currents (for this case +2 current). 
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Exercise 2. 


(3, 3) 


(-3, -3) 


Exercise 3. 


b) On the axis of the loop, 


where | is the current, b is the radius of the current loop, and x is 
the distance along the x-axis. You derive this result by using the 
Biot-Savart law analytically and noticing that, by symmetry, By and 
B, are zero. You can integrate the equation for dB,. 





(3, 3) 


c) y 


0.3 
0.6 


0.9 
1.2 
1.5 
1.8 


{(-3, -3) 


Bx 

6.28319 E-7 
6.74652 E-7 
8.86302 E-7 


2.46628 E-6 

-6.69063 E-7 
-1.78912 E-7 
—8.13398 E-8 
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Exercise 5. 


(3, 3) 





(-3, ~3) 


Exercise 7. 


a) Inside the field is zero (as well as this approximation will allow). 


Outside the field is the same as if it were from a single line current of 
magnitude 16 at the origin. 


>> 
b) The Biot-Savart law is $B*d& = kgl. Since any line integral inside the 


Cc 


‘wena 


cylinder will contain no current, the B field must be zero. Outside 
the cylinder, any line integral completely outside will contain all the 
current. Therefore, since the system is radially symmetric it will be 
the same magnitude and shape as the field from a single conductor at 

the origin. . 


The field inside both and outside both is zero. The field between the 
two is B = 2kg1/r where r is the radius from the origin. (This agrees 
with the Biot-Savart law and is similar to the problem of a Toroid.) 





Exercise 10. 


(7,7) 


(-7, -7) 


CHAPTER THREE 


Exercise 1. 


a) The trajectory moves away from the line charge but becomes 
straighter as the positron gets farther away (at constant z velocity). 


c) Note the angle is almost independent of distance but those really 
close get bent more. 


d) The positron is deflected toward the wire. 


e) Yes, but since the force goes as 1/r (instead of 1/r7), the orbits are 
(generally) more complicated than elliptical (satellite) orbits. 
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(40, 20) 


(-40, 0) 


Exercise 2. 


(+50, +5) 





(-10, -5) 





a) The electron will always be equally repelled by both line charges. 
Therefore, there are never any net forces on the electron and it will 
remain traveling straight along the z-axis at constant velocity. 


b 


— 


The electron will always be equally repelled by both line charges. 
Initially, the electron slows up because of a net force in the —x 
direction. If the electron did not have energy enough to pass the 
origin, it would turn around and go back the negative x-axis. Once 
the electron passes the origin, it will accelerate out the positive 
x-axis, 


~— 


c) Yes, and in such a way as to cause the x-axis crossing to be at the 


same point for all of the trajectories near the center. This can be 
thought of as a two-dimensional electron lens. You can also see the 
spherical aberration of trajectories near the line charges. 


Exercise 3. 
a) The positron circles both charges but the orbit does not close. 


b) The positron orbits just the charge at (3,0,0) but is perturbed by the 
field of the other charge. 


c) The positron orbits just the charge at (3,0,0) but is perturbed by the 
other charge so that its orbit changes. 


Exercise 5. 


a) The electron’s motion is unaffected by the field since the velocity is 
parallel to the magnetic field. 


b) The electron travels with a constant velocity in the z direction but 
spirals around an axis parallel to the z-axis. 


c) The electron just spirals around a circle in the x-y plane. 

d) The field collimates the electrons by forcing them to spiral around 
the magnetic field and not escape in some unwanted direction. 

Exercise 6. 

a) The positron when traveling into the converging field starts to travel 
in the y+z direction. Eventually, all of the positron’s velocity Is 
gone, and the positron is just traveling in the y-z plane. Then the 
positron starts coming back out. 


b) It decreases to zero, then changes direction. 


c) Energy is conserved, v2 stays constant. 
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Exercise 7. 


b) If the charge is positive, then the positron with V,>|E//|B] will move 
in the +y direction. When V,</EI/|BI, the positron will move in the 
-y direction. If the charge is an electron (a negative charge), the 
results are just reversed. 


“ae” 


c) A positively charged particle gets accelerated in the x direction (by 
> —> 

VyxBz). This acceleration then causes the velocity, Vy, to grow 
which, in turn, creates an increase in Vy: If the charge is negative, 


Vy decreases and V\, increases. 


d) The z motion is unaffected. The xy motion is just that already 
discussed. 


e) If there is no change in velocity of the particle, the forces due to E 
and B must be equal. 


> >> 

gE = qVxB 

Ey = V,B, 

Vy = Ey/Bz with no requirements on V, 
- Exercise 9. 

-> 

F=ma 

> > 

a= F/m 

> -> re 

V = Viinitial) * a 

> > = >.9 

S= S(initial) + V(initiat) t+ 1/2 at 


since a= ay = Eq/m 
Sy ~ Sy (initial) = Vy (initial) t + (1/2)(Eq/m) t? 
Sy = Sy (initial) + Vx(initial) t 
Sz = Sz(initial) + Vz(initial) t 


For the example in the text: 


The initial conditions are ee Sy = S$, = 0 and initial velocities are Vy = 
V, =Qand Vy =4 

Sy = 1/2 Eq/mt? =1/20? 

Sy = 4t 


S, =0 
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Exercise 10. 


) _mVv__v. 
a r= ciel B 


—_> 
when |B| is measured in units of q/m. 


mp 


b) T 


li 
— fo. 
Il 
NO 
po 
ONS 


= 


in normalized units. 


CHAPTER FOUR 


Exercise 1. , 


a) The surface integral equals +1 times the fength of the box along the 
z-axis. The line integral is zero. 


b) The surface integral equals +2 times the length of the box along the 
z-axis. The line integral is zero. 


c&d) Both the surface and line integrals are zero. 


Exercise 2. 


a) The surface integral is zero. The line integral is +1 and independent 
of the contour. 


b) The surface integral is zero. The line integral is +2 and independent 
of the contour. 


c&d) Both the surface and line integrals are zero. 


Exercise 3. 


a) The surface integral is +3 times the length of the box along the 
z-axis. The line integral is zero. 


b) The surface integral is +3 times the length of the box along the 
z-axis. The line integral is 2. 


c&d) Both the surface and line integrals are zero. 


Exercise 4. 


a) 8. The surface integral is +1 times the total area of the box perpen- 
dicular to the z-axis. 


b) It is caused by a plane of charge along the x-y axis. 


c) zero 
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Exercise 5. 
a} zero 
c) 0,4,0 


d) it is caused by a plane (defined by the y-z plane) of current traveling 
in the y direction. 


Exercise 6. 
b) zero 


> -> A => 
c) $F°d& = $rf(r)-d& 
25 8 
If we pick any circle around the origin, then 7 1 d&. Therefore the 


integral is zero. 


Exercise 7. 
a) zero 


b) There are no sources or sinks for the fluid, that is, there are no 
places where fluid is being created or destroyed. 


c) zero 


d) -.75 


CHAPTER FIVE 


Exercise 1. 


a) 


# OF ITERATIONS = 60 
0 6 0. 0 0 0 0 9 0 0 0 


2 2oF 165 Bn & FoH 165 10% Bt Be6 363 4 

4 304 209 2e% 206 208 Fod 4 Sol 624 8 

6 § Meo2 JeoF BoF BoP 406 S$oF Fe® DeA 12 
8 603 Set 465 463 466 5-64 608 809 32109 16 
100—Cti«? SoS 408 @6 499 So®? Fok 9e6 13e3 20 
8 6e3 Sel 205 403 06 SoH 608 8.9 11-9 16 
6 5 Me2 SeF 3oF 309 266 FoF Fe eA 12 
4 3eoM 209 B2eF B26 2-68 363 & S.-i 604 §& 


2 $oF 10985 oS 104 105 eF Beh BH BoD 4 
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c) The pattern is symmetric across the x-axis (i.e., across the line 
joining the +10 and +20 potential points.) 


Exercise 2. 


a) 
¢ OF ITERATIONS = 186 

0 0 0 0 0 0 0 0 0 0 0 

2 let OFS Get O Oo! O63 OfF8 166 266 & 

4 Ze 103 Oeo6 Oo Oo OoF 369 3% Set 8 

6 3e? 2el tet Go? GeS 366 3 Se-1 6el 12 
8 409 2.6 128 1! n2 Be1 368 GeS 1008 16 
10 $e6 Jot %$eF bel Fed Be3 Bel Fe® 1349 20 
8 4.9 2.68 1-5 1! 10-2 203 3c8 6653 10-8 56 
6 307? Be! tol Ge? O98 6 3 $e-1 68-3) 12 


4 Seh $e3 OF6 Bod Ooh OF $1069 3ed Set 8 


2 let OeS Oot 0 Ov} Oc3 Oo8 86 266 4 


b) The pattern is symmetric across the x-axis (i.e., across the line 
joining the +10 and +20 potential points.) 


Exercise 5. 


a) The symmetries are along the x- and y-axes and also along the two 
diagonals of the square. 


b) 

¢ OF ITERATIONS = 190 
0 
1 108 


2 2-7 Jed 


3 30S fet 406 


4 4ce4 46558 5022 326 


§ $e2 SA 5.7 $e9 6.2 


t 608 60% 626 6e7F nF 608 629 


8 4e6 72-2 =F 6.9 6099 6.9 7 3 


10 SoS BPeS BPe@ FeO Tel F 3 Fei Pel Bel 
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Exercise 6. 


a) The symmetries are along the x- and y-axes and also along the two 
diagonals of the square. 


b) 
¢ OF ITERATIONS © S2 
0 


3 308 403 $22 


6 604 668 Fe2® 2e6 8&8 Se3 

7 Fe2 Ted TeT? 8 Se3 Seb Be8 

8 8 8 6-2 Se4 8-6 6-6 9 98 

9 Se6 Goh 808 8-6 BeS 9D 9e1 Ded FoS 

i00— («9 SoG 8x6 8.2% Se8 9 Ge2 Fed Fe Bee 


Exercise 7. 


a) The symmetries are along the x- and y-axes and also along the two 
diagonals of the square. 


(9,9) 

0 

2 0 

9 5 0 

21 14 7 0 

37. 27 #18 9 0 

60 46 34 22 = «11 0 

100 73 55 39 26 —= 13 0 

100 100 78 59 43 28 ~— 14 0 
100 100 100 80 62 45 29 14 O- 

100 100 100 100 ~~ 81 62 46 30 15 0 
(0,0) (9,0) 


Exercise 8. 
a) The pattern is symmetric across the x-axis. 
b) The potential along the y-axis is zero. 
c) Since we know the potential along the y-axis is zero, we can use the 


axis as one of the boundaries to the problem. In this way, the region 
to be solved only contains one charge. 





